Kepler-47: A Transiting Circumbinary Multi-Planet System 
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We report the detection of Kepler-47, a system consisting of two planets orbit- 
ing around an eclipsing pair of stars. The inner and outer planets have radii 
3.0 and 4.6 times that of the Earth, respectively. The binary star consists of a 
Sun-like star and a companion roughly one-third its size, orbiting each other 
every 7.45 days. With an orbital period of 49.5 days, eighteen transits of the 
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inner planet have been observed, allowing a detailed characterization of its or- 
bit and those of the stars. The outer planet's orbital period is 303.2 days, and 
although the planet is not Earth-like, it resides within the classical "habitable 
zone", where liquid water could exist on an Earth-like planet. With its two 
known planets, Kepler-47 establishes that close binary stars can host complete 
planetary systems. 

The extremely precise and nearly continuous observations provided by the Kepler spacecraft 
(1) has enabled the detection of over 2300 planet candidates (2, 3), and over 2100 eclipsing 
binary stars (4, 5). A synergy of these efforts has helped establish the class of circumbinary 
planets, which are planets that orbit around a pair of stars (6-8). Their detection has led to a 
revitalized effort to understand planet formation around binary stars (9, 10). A circumbinary 
planet can reveal itself in two ways. If the planet's orbital plane is favorably aligned, the planet 
may transit across one or both of the stars, causing a small decrease in the amount of light from 
the system. If the planet is sufficiently massive and close, the planet can perturb the stellar 
orbits (11). The most readily observable manifestation of this perturbation is a change in the 
times when the eclipses occur. 

In contrast to a single planet orbiting a single star, a planet in a circumbinary system must 
transit a "moving target." As a consequence, the time intervals between the transits as well as 
their duration can vary substantially. The transits can deviate from having a constant period by 
up to several days and can vary in duration by several hours. These transit timing and duration 
variations can be taken to be the signature of a circumbinary planet because no other known 
mechanism can cause such effects. Modeling of the timing and duration changes can be used to 
precisely determine the orbits of the planet and stars (6-8). 

Kepler observations of the binary star system Kepler-47 (KIC 10020423, also KOI-3154) 
show primary eclipses (smaller and fainter star blocking the brighter, more massive "primary" 
star) every 7.45 days with a depth of 13%. Also present are secondary eclipses with a depth of 
0.8% (Fig. 1) and a quasi-periodic modulation in the out-of-eclipse regions of ~ 2 — 4% caused 
by star-spots on the primary star (12). The Kepler data span 1050.5 days and visual inspections 
of the light curve revealed the signals of two candidate circumbinary planets, with periods of 
~ 50 and ~ 303 days. Three transits of the longer-period candidate (hereafter the outer planet) 
were readily apparent, but those of the shorter-period candidate (hereafter the inner planet) were 
more difficult to find because of their shallower depth. Using the predictions of a preliminary 
model of the system as a guide, a total of 18 transits of the inner planet were detected. The 
transits have timings that can deviate from strict periodicity by up to several hours and their 
durations vary significantly, strongly suggesting their origins are from circumbinary planets 
(Fig. 2). All of these events are transits over the primary star. 

To characterize the stellar orbit, we obtained Doppler spectroscopy of the system (12, Fig. 
1). Radial velocity variations of the primary star were readily detected, but the secondary star 
is too faint to have been measured. Usually when the radial velocity measurements of only one 
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component in a spectroscopic binary are available, the masses of the stars cannot be uniquely 
determined. However, the transit times and durations provide constraints on the geometric 
configuration of the stellar orbits and specify the stellar mass ratio, which in combination with 
the primary's radial velocities, allow both stellar masses and the physical scale to be determined. 

To determine the system parameters, we used a photometric-dynamical model (13) similar 
to that used for the four previously known transiting circumbinary planets (6-8). This model 
assumes spherical bodies interacting via Newtonian gravity (12), and is used to fit the radial 
velocity data and the Kepler time- series photometry. We determined the stellar masses as de- 
scribed above, and the relative sizes of the bodies from the eclipses and transits in the light 
curve. Information on the inclination, eccentricity, and mutual inclination of the planetary or- 
bits is also implicit in the combination of photometric and radial-velocity data. Gravitational 
perturbations caused by the planets on the stars and on each other could, in principle, also con- 
strain the masses, but for Kepler-47 the expected masses of the planets are too small to create 
a measurable effect over the time span of our data. The small radii of the transiting objects 
strongly suggests they are of planetary mass (Table 1); dynamical considerations described be- 
low make this conclusion secure. 

The inner planet, Kepler-47 b, is the smallest transiting circumbinary planet yet detected, 
with a radius of 3.0±0. 1 Earth radii. Its mass is too small to be directly measured, but a 3a upper 
limit of 2 Jupiter masses has been determined based on the nondetection of timing variations 
of the stellar orbit (12). Because the planet's mass is unknown, its density is also unknown 
and it is not possible to distinguish between a rocky composition and a more volatile-enriched 
composition. We can make a plausible mass etimate by using both an empirical mass-radius 
relation based on transiting exoplanets (14) and a limited empirical mass-radius relation for 
planets in the Solar System (15) yielding ~ 7 — 10 Earth masses or ~ 0.4 — 0.6 Neptune 
masses. The planet's 49.5-day orbital period is 6.6 times the period of the stellar binary. This 
is ~ 77% longer than the critical period (28 d) within which the planet would be susceptible to 
dynamical instability due to interactions with the stars (16). While this 77% margin is notably 
larger than for the other known transiting circumbinary planets, i.e., 14%, 21%, 24%, 42% for 
Kepler 16, 34, 35, and 38, respectively, the planet is still somewhat close to the instability limit, 
a feature shared by all known transiting circumbinary planet systems. 

The outer planet, Kepler-47 c, has a radius of 4.6 ± 0.2 Earth radii, making it slightly larger 
than the planet Uranus. As before, the planet's mass is too small to be measured directly, and 
we derived a 3a upper limit of 28 Jupiter masses (12). Based on its radius, we find a plausible 
mass of ~ 16 — 23 Earth masses or ~ 0.9 — 1.4 Neptune masses, using these empirical mass- 
radius relations (14, 15). With only 3 transits currently available, the outer planet's orbital 
eccentricity is poorly constrained. A perfectly circular orbit would fit the data, and a low- 
eccentricity orbit seems plausible given the low eccentricity of the stellar binary (e = 0.023) 
and of planet b (e < 0.035). The photometric-dynamical model provides only an upper limit on 
the eccentricity, e < 0.4 with 95% confidence, and the requirement of long-term stability only 
rules out eccentricities larger than 0.6 (12). 

Due to the orbital motion of the stars, the outer planet is subject to variations in the incident 
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stellar flux (i.e. insolation), even if the planet's orbit is circular (Fig. 3). The average insolation 
is similar to the amount the Earth receives from the Sun: for a circular orbit it is 87.5% of the 
Sun-Earth insolation, and varies by ~ 9%. This places Kepler-47 c well within the classical 
"habitable zone", defined as the range of distances from the host star(s) where liquid water 
could persist on the surface of an Earth-like planet (17). While Kepler-47 c is probably a gas 
giant and thus not suitable for life, its location is notable as it demonstrates that circumbinary 
planets can exist in habitable zones. Although the definition of the habitable zone assumes a 
terrestrial planet atmosphere which does not apply for Kepler-47 c, large moons, if present, 
would be interesting worlds to investigate. 

A 0.2% deep transit-like event is present at time 2,455,977.363 (BJD) that is not caused by 
either of the two planets. A search for additional transits has revealed several more tentative 
transit events (72), but we caution that the star is faint (the Kepler magnitude is 15.178), there 
are large modulations due to star-spots, and the data contain correlated "red" noise, making 
small, non-periodic transit detection challenging. The marginal evidence at the present time is 
insufficient to place confidence on any additional candidate planet(s). 

The primary star is similar to the Sun in both mass and radius, and dominates the luminosity 
of the binary system, having 60 times the bolometric luminosity of the secondary star (or 176 
times the brightness in the Kepler bandpass). A spectroscopic analysis gives an effective tem- 
perature of 5640 ± 100 K for the primary star (Table S2), with a metallicity slightly less than 
solar ([M/H]= —0.25 ± 0.08 dex). The star's rotation period as determined from the star-spot 
modulation in the light curve (12) is only 4% longer than the orbital period, suggesting that the 
spin and orbital angular momenta have been synchronized by tidal interactions. Supporting this 
interpretation, the obliquity of the primary star (the angle between the spin and orbital axes) 
must be smaller than about 20°, based on the observable effects of the secondary star eclips- 
ing star-spots on the primary star (12, 18-21). Star-spot crossings also perturb the shape and 
depth of the primary eclipses, leading to systematic trends in the eclipse times, and limit the 
precision with which one can infer the planets' masses. In addition, the loss of light due to star- 
spots causes eclipses to appear slightly deeper than they would for an unspotted star, biasing 
the determination of the stellar and planetary radii too high by a few percent. 

With Kepler-47 b and c, there are six confirmed transiting circumbinary planets currently 
known. Their orbital periods relative to their host binary stars show no tendency to be in res- 
onance, and their radii are Saturn-size and smaller. Given that Jupiter-size planets are easier 
to detect, their absence in the Kepler data suggests that the formation and migration history of 
circumbinary planets may disfavor Jupiter-mass planets orbiting close to the stars, in accord 
with (22). 

The planets in Kepler-47 are expected to have formed much farther out than their present 
orbits, at locations where the conditions for the formation of giant planets are more favor- 
able (9, 10). The planets have likely migrated to their current orbits as a result of interactions 
with the circumbinary disk. The multiplicity and coplanarity of the orbits strengthens the argu- 
ment for a single-disk formation and a migration scenario for circumbinary planetary systems. 
However, unlike orbits around a single star, the environment around a binary star is much more 
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dynamic and tends to augment planet-planet interactions. The relatively large distance between 
the orbits of the inner and outer planets in the Kepler-47 system is consistent with requirements 
for dynamical stability (23). 

The previously detected transiting circumbinary planet systems show no evidence for more 
than a single planet. The multi-planet nature of the Kepler-47 system establishes that despite 
the chaotic environment around binary stars, planetary systems can form and persist close to the 
binary, and invites a broader investigation into how circumbinary planets compare to planets and 
planetary systems around single stars. 
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parameter 


best fit 


la uncertainty 


Bulk Properties 


Mass of Star A, M A (M ) 


1.043 


0.055 


Mass of Star B, M B (M ) 


0.362 


0.013 


Radius of Star A, R A (Rq) 


0.964 


0.017 


Radius of Star B, R B (R e ) 


0.3506 


0.0063 


Temperature of Star A, T eff ^ (K) 


5636 


100 


Temperature of Star B, T cff B (K) 


3357 


100 


Luminosity of Star A, L A (L Q ) 


0.840 


0.067 


Luminosity of Star B, L B (L e ) 


0.014 


0.002 


Radius of Planet b, R b (-R©) 


2.98 


0.12 


Radius of Planet c, R c (R e ) 


4.61 


0.20 


Stellar Orbit 


Semimajor Axis, a A B (AU) 


0.0836 


0.0014 


Orbital Period, P AB (day) 7.44837695 


0.00000021 


Eccentricity, e AB 


0.0234 


0.0010 


Argument of Periapse, uo AB (Degrees) 


212.3 


4.4 


Orbital Inclination, % x (deg) 


89.34 


0.12 


Planet b Orbit 


Semimajor Axis, a b (AU) 


0.2956 


0.0047 


Orbital Period, P b (day) 


49.514 


0.040 


Eccentricity (95% conf.), e b 


< 0.035 




Orbital Inclination, i b (deg) 


89.59 


0.50 


Mutual Orbital Inclination, I b (deg) 


0.27 


0.24 


Planet c Orbit 


Semimajor Axis, a c (AU) 


0.989 


0.016 


Orbital Period, P c (day) 


303.158 


0.072 


Eccentricity (95% conf.), e c 


< 0.411 




Orbital Inclination, i c (deg) 


89.826 


0.010 


Mutual Orbital Inclination, I c (deg) 


1.16 


0.46 



Table 1 : A summary of the results for the photometric-dynamical model. For brevity some 
of the fitting parameters are not listed here. See Table S5 for a complete listing of fitting param- 
eters. 
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Fig. 1: Light curves and velocity curve data with model fits. Top: Normalized and detrended 
flux is plotted versus orbital phase for the primary and secondary eclipses, along with the binary 
star model. Middle: The radial velocities of the primary star and the best-fitting model are 
plotted versus the orbital phase. The expected radial velocity curve of the secondary star is 
shown with the dashed line. Bottom: The normalized and detrended flux near five representative 
transits of the inner planet and all three transits of the outer planet are shown. See Figs. S13, 
SI 4, and S15 for plots of all 18 transits of the inner planet and plots of the residuals of the 
various model fits. 
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Fig. 2: Planetary transit time and duration variations. Left: The observed minus expected 
times of transit computed from a linear ephemeris are shown versus time (an "O-C" curve). 
The triangles show the measured deviations, and the filled circles are the predictions from the 
photometric-dynamical model. Four transits of the inner planet occurred in data gaps or regions 
of corrupted data. Top right: The O-C values of the inner planet are shown as a function 
of the binary phase, where the primary eclipse occurs at phase 0.0 and the secondary eclipse 
is at phase 0.487. Two cycles have been shown for clarity. The solid curve is the predicted 
deviation assuming a circular, edge-on orbit for the planet. The lateral displacement of the 
primary near the eclipse phases is minimal and therefore the deviation of the transit time from a 
linear ephemeris is near zero. The primary is maximally displaced near the quadrature phases, 
so transits near those phases show the most offset in time. Bottom right: The durations of the 
transits for the inner planet (filled circles) and the outer planet (filled squares) as a function of 
the orbital phase of the binary. The solid curves are the predicted durations assuming a circular, 
edge-on orbit for the planet. At phases near the primary eclipse, the planet and the primary star 
are moving in opposite directions, resulting in a narrower transit. At phases near the secondary 
eclipse, the planet and the primary star are moving in the same direction, resulting in a longer 
transit. The outer planet is moving slower than the inner planet, resulting in longer transits at 
the same binary phase. 
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Fig. 3: The time- varying insolation S received by Kepler-47 c, for different assumed eccen- 
tricities. The insolation is in units of the Solar luminosity at a distance of 1 AU (Ss un = 1368 W 
m~ 2 ). The upper panels are for a zero eccentricity orbit and highlight the insolation variations 
caused by the 7.4-day orbit of the binary. The lower panels show eccentricities of e = 0.0, 0.1, 
0.2, and 0.4 (colored black, green, blue, and red, respectively), and illustrate the longer time- 
scale variations. The dotted lines mark the limits for the inner and outer edges of the habitable 
zone, following the prescription in (24) for the onset of a runaway greenhouse effect and the 
maximum greenhouse effect. 
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Supporting Online Material (SOM) 



We provide additional details regarding the detection and characterization of Kepler-47 
in this supplement. §l.ll gives alternate designations and other information for Kepler-47. 
§1.21 discusses the Kepler data preparation and detrending. §1.31 discusses how the rota- 
tional period of the primary star is derived. §1.4l discusses the ground-based spectroscopic 
observations. §1.51 describes how the effective temperature, gravity, and metallicity of the 
primary were measured. §1.61 gives an overview of how the times of mid-eclipse for the 
primary and secondary eclipses were measured. §1.71 discusses the effects of star-spots on 
the measurement of the eclipse times and other parameters. §1.81 presents measurements 
of the transit times and the detection of a transit event possibly due to a third planet. §1.91 
gives a full discussion of the photometric-dynamical model. §1.101 presents a discussion 
of independent light curve modeling done with the ELC code. §1.111 discusses how upper 
limits on the masses of the planets were derived. §1.12l considers the long-term stability of 
the planetary orbits. §1.131 gives a comparison of the stellar properties with evolutionary 
models. §1.14l presents details of the habitable zone in Kepler-47. 

1 Materials and Methods 

1.1 Alternate designations, celestial coordinates, and apparent magnitudes 

Kepler-47 appears in the Kepler Input Catalog (25, KIC) as KIC 10020423. Other designations 
include Kepler Object of Interest KOI-3154 and 2MASS J1941 1149+4655136. The J2000 ce- 
lestial coordinates given in the KIC are a = 19 h 41 m llf 501, 5 = +46°55'13"69, and the appar- 
ent magnitudes are r = 15.126 and Kp = 15.178. 

1.2 Kepler data preparation and detrending 

In this study we make use of data from Kepler Quarters Ql through Q12 (May, 2009 through 
late March, 2012). We used the "simple aperture photometry" (SAP) provided by the Kepler 
pipeline and available at the Mikulski Archive for Space Telescopes (MAST). The Kepler SAP 
light curves show instrumental trends (2(5), so further processing is necessary. The detrending 
must be done for each Quarter separately since the object appears on a different detector module. 
The amount of detrending needed depends on the specific task. When modeling the eclipses and 
transits, a fairly aggressive detrending is used where both the instrumental trends and the spot 
modulations are removed. In this case, the eclipses and transits are masked out, and a high order 
cubic spline is fit to short segments whose end points are usually defined by gaps in the data 
collection due to monthly data downloads, rolls between Quarters, or spacecraft safe modes. 
The segments are normalized to the spline fits, and the segments are reassembled. The SAP 
light curves and the detrended light curve with the spot modulation removed are shown in Fig. 
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|S4| Other tasks such as spot modeling require much less aggressive detrending, in which case 
low-order polynomials are used to stitch together different segments across the Quarters. 

The time difference between the last Q12 observation and the first Ql observation is 1050.51 
days. During that interval, Kepler was collecting data 92.55% of the time, and 44389 cadences 
out of the 47580 in total were flagged as good (SAP_QUALITY=0), for a duty cycle of 86.34%. 
Not all observations with SAP_QUALITY>0 are necessarily useless, depending on the purpose, 
so the 86.34% duty cycle is a lower limit. 

1.3 Rotational period from star-spot induced stellar variations 

Fig. [S5] shows closer-in views of the light curves from Ql, Q5, and Q9. A modulation of up to 
3% in the out-of-eclipse regions due to star spots rotating into and out of view is evident. This 
modulation has a period that is close to, but not exactly equal to the eclipse (e.g. orbital) pe- 
riod. Fig. [S6] depicts the autocorrelation of the cleaned detrended light curve, after the primary 
and secondary eclipses were removed and replaced by the value of the mean light curve with a 
typical random noise. The autocorrelation reveals clear modulation with a period of about 7.8 
days. Presumably, the clock behind the modulation is the stellar rotation of the primary, which 
has brightness variation due to inhomogeneous distribution of stellar spots (as the primary star 
dominates the light in the Kepler bandpass, we assume it is the source of the modulation). To 
obtain a more precise value of the stellar rotation we measured the lags of the first 12 peaks of 
the autocorrelation and fitted them with a straight line as shown in Fig.[S7J From the slope of the 
fitted line we derived a value of 7.775 ± 0.022 days as our best value for the stellar rotation pe- 
riod. This period is slightly longer than the orbital period of 7.448 days. It is interesting to note 
that the transition between synchronized and unsynchronized binaries for pre-main sequence 
and young stars appears between 7 and 8 days, as depicted by (27). 

1.4 Spectroscopic observations 

We observed Kepler-47 four times with the High-Resolution-Spectrograph (28, HRS) at the 
Hobby-Eberly Telescope (HET). Spectra with a resolving power of R = 30, 000 were obtained 
on UT 2012, April 23, May 18 & 20 and June 5. We used the "600g5822" setting of HRS that 
delivers a spectrum from 4814 to 6793 A. The data were reduced with our own HRS reduction 
script using standard IRAF routines. We selected a total exposure of 3600 seconds per spectrum 
(divided into three sub-exposures of 1200 seconds each to facilitate cosmic-ray removal). The 
signal-to-noise (S/N) levels of the HRS spectra range from 30:1 to 55:1 at 5500 A, depending 
on seeing conditions. Adjacent to every visit to Kepler-47 we also observed the Kepler field 
standard star HD 182488 to be used for the radial velocity determination. 

In addition to the HET observations, we observed Kepler-47 six times using the Tull Coude 
spectrograph (29) at the Harlan J. Smith 2.7m telescope (HJST). The data were obtained with 
our standard instrumental setup that covers the a wavelength range of 3760-10,200 A and uses 
a 1.2 arcsecond slit that yields a resolving power of R = 60, 000. We obtained data during the 
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nights of UT 2012, May 1, 2, 4-6 and on June 26. Exposure times ranged from 3600 to 4800 
seconds (again divided in 1200 second sub-exposures) and the S/N is typically around 14:1 at 
5500A. Each of these nights we also observed HD 182488 to serve as a radial velocity standard. 
The data were reduced with our own reduction scripts using standard IRAF routines. After 
some experimentation, it was discovered that better measurements of the radial velocities were 
obtained from spectra that did not have the sky background subtracted. 

An additional spectrum of Kepler-47 was obtained using the 10 m Keck 1 telescope and the 
HIRES spectrograph (30). The spectra were collected using the standard planet search setup 
and reduction (31). The resolving power is R = 60, 000 at 5500 A. Sky subtraction, using the 
"C2 decker" was implemented with a slit that projects to 0.87 x 14.0 arcsec on the sky. The 
wavelength calibrations were made using Thorium- Argon lamp spectra. 

The radial velocities of Kepler-47 were measured using the "broadening function" technique 
(32). The broadening functions (BFs) are rotational broadening kernels, where the centroid of 
the peak yields the Doppler shift and where the width of the peak is a measure of the rotational 
broadening. The BF analysis is often better suited for measuring radial velocities of binary stars 
in cases where the velocity difference between the two stars is small compared to the spectral 
resolution. A high quality spectrum of a slowly rotating star is needed for the BF analysis, 
and for this purpose we used observations of HD 182488 (spectral type G8V) taken with each 
telescope+instrument combination. The derived radial velocities are insensitive to the precise 
spectral type of the template, as similar radial velocities are found when using templates of early 
G to late K. The adopted template radial velocity was —21.508 km s _1 (33). 

We prepared the spectra for the BF analysis by normalizing each echelle order to the local 
continuum using cubic splines, trimming the low signal-to-noise ends of each order, and merg- 
ing the orders by interpolating to a log-linear wavelength scale. The wavelength ranges used 
for the final BF analysis was 4830-5770 A for the HET spectra and 5138-5509 A for the HJST 
spectra. Fig.[S8]shows example BFs from HET and HJST spectra. The spectrum is single-lined, 
as only one peak is evident in the BFs from the HET. Some simple simulations were performed, 
and non-detection of a second star in the HET spectra indicates the secondary star is >10 times 
fainter than the primary star, consistent with the expectations based on the eclipse depths, where 
a flux ratio of ~ 1/176 is expected. In the case of the HJST, two peaks are apparent. However, 
one of the peaks is due to the sky background since it is stationary in velocity, and changes 
strength relative to the other. The FWHM of the BF peaks were consistent with the instrumen- 
tal broadening, which indicates the rotational velocity of the primary is at best only marginally 
resolved. 

Gaussian functions were fit to the BF peaks to determine the relative Doppler shifts. The 
appropriate barycentric velocity corrections were applied and the contribution of the template 
radial velocity was removed, thereby placing the radial velocities on the standard IAU radial 
velocity scale defined by (33) and (34). The Keck HIRES pipeline automatically produces radial 
measurements for single stars on the IAU scale, accurate to 0.1 km s _1 or better (34). Having 
established Kepler-47 as a single-lined binary, we simply adopted the pipeline measurement. 
The radial velocity measurements for all 11 observations are given in Table SI. 
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1.5 Spectroscopic parameters 

The effective temperature Teg, surface gravity \ogg, the metallicity [m/H], and the rotational 
velocity V mt sin i of the primary were measured using the Stellar Parameter Classification (SPC) 
code (35). SPC uses a cross -correlation analysis against a large grid of model spectra in the 
wavelength region 5050 to 5360 A. Since all of the absorption lines in this region are used, the 
SPC analysis is ideal for spectra with low signal-to-noise. The first three HET observations were 
combined to yield a spectrum with a signal-to-noise ratio of m 116 in the order containing the 
Mg b features near 5169 A (the fourth HET observation had relatively high sky contamination 
and was not used). The derived spectroscopic parameters are given in Table S2. 

1.6 Stellar eclipse times and corrections 

The times of mid-eclipse for the primary and secondary eclipses in Kepler-47 were measured 
using the technique described in (7). For completeness we give most of the details here as well. 
Given an initial linear ephemeris and an initial estimate of the eclipse widths, the data near the 
eclipses were isolated and locally detrended using a cubic polynomial with the eclipses masked 
out. The detrended data were then folded on the linear ephemeris and an eclipse template was 
made by fitting a cubic Hermite spline. The Piecewise Cubic Hermite Spline (PCHS) model 
template was then iteratively cross-correlated with each individual eclipse event to produce a 
measurement of the time at mid-eclipse. After each iteration, a new PCHS model was produced 
by using the latest measured times to fold the data. Fig. [S9] shows the folded eclipse profiles and 
the final PCHS models. The fits are generally good, although there is increased scatter near the 
middle of the primary eclipse, presumably due to the effects of star spots. Table S3 gives the 
eclipse times. The cycle numbers for the secondary eclipse are not exactly half integers because 
the orbit is eccentric. 

The eclipse times were fitted with a linear ephemeris and the Observed minus Computed (O- 
C) residual times were calculated and are shown in Fig. lSlOl For the primary, there are coherent 
deviations of up to two minutes. While not strictly periodic, there is a quasiperiod of « 178 
days seen in a periodogram (Fig. IS 1 ll) . This modulation is most likely a beat frequency between 
the stellar rotation and the binary motion, similar to what is observed for Kepler- 17 (20). If 
the secondary passes in front of a big spot during the primary eclipse, the spot anomaly will 
introduce a shift on the eclipse timing since the projected stellar disk of the primary on the 
sky will no longer have a symmetric surface brightness distribution. The shift of the eclipse 
time will depend on the size and position of the spot and the position on the eclipse chord. A 
long-lasting spot will introduce shifts in consecutive eclipses, but the shift will change with 
time since the spot will be at a different position on the eclipse chord at each eclipse. More 
specifically for this system, a spot with a period of rotation of 7.775 days will effectively recede 
on the transit chord 360°(7.4484 d - 7.775 d)/7.4484 d = -15.79° each eclipse. In order to 
come back to the exact same position, and hence complete a full cycle in the O-C diagram, 
(360°/15.79°)P O rb = 22.8P orb = 170 days will be needed, which is close to the period of the 
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observed signal. In reality, the spots also change with time and may also drift in latitude over 
time, so the signal near the beat frequency is blurred somewhat. 

There is a correlation between the O-C residual time of the primary eclipse and the local 
slope of the out-of-eclipse portions of the SAP light curve during the eclipse, as shown in Fig. 
IS 121 A large negative slope in the light curve surrounding an eclipse indicates a dark spot is 
rotating into view. The "center of light" of the primary will be shifted to the opposite side of 
the stellar disk, resulting in a slightly later time of mid-eclipse. Likewise, a large positive slope 
surrounding an eclipse indicates a dark spot is rotating out of view, which results in a slightly 
earlier time of mid-eclipse. Finally, when the slope is near zero, the spots are centered on the 
stellar disk, and no change in the eclipse time is seen. A linear function was fitted to the data 
in Fig. IS 121 and the times of primary eclipse were corrected. The O-C diagram resulting from 
these corrected times (Fig. lSlOl) has much less scatter. No periodicities are evident (Fig. IS 1 ll) . 

The best-fitting ephemerides for the corrected primary eclipse times and the secondary 
eclipse times are 

P A = 7.44837605 ± 0.00000050 d Kepler-47 primary 

P B = 7.44838227 ± 0.00000342 d Kepler-47 secondary 

T (A) = 2, 454, 963.24539 ±0.000041 Kepler-47 primary 

T (B) = 2, 454, 959.426986 ±0.000277 Kepler-47 secondary 

The difference between the primary and secondary periods is 0.52 ± 0.30 seconds, with the 
secondary period being longer. 

1.7 The effect of star-spots on the eclipses: possible biases and spin-orbit 
alignment. 

Star-spots cause the light curve to exhibit modulations that can be used to measure the rotation 
period of the primary star. Star-spots can also affect the determination of certain system param- 
eters. It has been shown that there is a correlation between the eclipse timing variations and 
the local slope of the stellar flux variations at the times of the eclipses. In order to confirm that 
star-spots are the main cause of the eclipse timing variations, and to evaluate their impact on 
our ability to measure the size of the secondary star, we attempt to model the effect of spots on 
individual eclipses (21, 36, 37). 

The data from each primary eclipse are isolated by keeping only 3 hours of observations 
before and after the eclipse. The out-of-eclipse part of each dataset is then fitted with a linear 
function. The fit is subtracted from the data, then the data are normalized so the out-of-eclipse 
flux is equal to unity. The detrended eclipse light curves are folded with a linear ephemeris, and 
this folded light curve is fitted with a standard model for light loss due to a dark body passing in 
front of a limb darkened star (38). This no- spot model has only four free parameters: squared 
radius ratio (R b /Ra) 2 , impact parameter b, normalized semimajor axis for the secondary orbit 
Ra/o-b, and a linear limb darkening coefficient u\. 
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The effect spots have on individual primary eclipses is manifest in two ways: the depth of 
each eclipse changes since it is measured relative to the changing stellar flux, and the shape 
of each eclipse is distorted which leads to a shift in the measured mid-eclipse time. Visual 
inspection of the eclipse residuals shows that this last effect can be well-modeled in most cases 
by adding just one large star-spot on the surface of the primary star. Since the rotation of the 
star happens on a longer time-scale than the eclipse itself, we held the position of the star-spot 
fixed during each individual eclipse. The latitude of the spots cannot be well constrained with 
single eclipse events, so we fix the position of the spot so that the center of the secondary star 
trajectory intersects the center of the spot. Our spot model adds five additional parameters: 
three parameters describe the spot itself - its angular radius, the flux contrast (related to the 
spot temperature), and the position along the eclipse chord. The fourth parameter is the out-of- 
eclipse flux, which corrects for the depth variations. Finally, the time of mid-eclipse is free. We 
set up a pixilated model of the star with a circular spot, in which the flux is calculated as the 
surface integral of the intensity of the visible hemisphere of the star. 

The best-fitting model for several consecutive eclipses is compared with a no-spot model in 
Fig. IS 1 3 L showing how the model captures the essential effect of spots on the eclipses. For every 
eclipse we obtained a new eclipse time, and fitted a linear ephemeris to these times. The scatter 
was found to be substantially reduced from the initial timings (see Fig. IS 101 upper panel), by 
a factor of 30%, which shows that indeed the scatter is due to spots. The improvement on the 
scatter is similar to the one obtained through the local slope correction, so this serves as a good 
consistency check. 

Our model also estimates the fraction of the star covered by spots at the time of each eclipse. 
This quantity is not very precise, but can help us estimate the effect of spots on the measurement 
of the eclipse depth and hence the radius ratio Ra/Rb- We divide the square of the radius ratio 
from the spot model by each observed local out-of-eclipse flux to mimic the apparent depth that 
one would obtain by fitting each eclipse individually. The results are plotted in Fig. IS 141 where 
one can clearly see how the depth of each eclipse changes with time. The variations seem to 
have a time-scale similar to the uncorrected eclipse timing variations, which is expected since 
the scatter in both are due to spots. A variation with the observing season is also apparent, 
which is a clear indication that there are different levels of Quarterly contamination. With these 
eclipse depths we can estimate the inferred secondary star radius Rb from each eclipse, using a 
fixed R A from Table 1 (see Fig. IS 141) . The values obtained do not have a large scatter, and they 
all agree within la with the value obtained from the photometric-dynamical model fit. Thus 
the correction to the secondary star radius because of the presence of spots is not significant, 
although a slightly smaller radius is favored. 

We can also use these spot models to gain information about the obliquity of the sys- 
tem (18-21, 39). In Fig. IS 1 31 we can clearly see how the spot model shows that a spot is 
moving backwards with each eclipse, which means that the spot trajectory is contained within 
the boundaries of the eclipse chord. This backwards movement makes it seem as if the star is 
rotating backwards (retrograde) very slowly, but this is simply a stroboscopic alias effect. The 
spot appears to move backwards because the star's rotation period is slightly longer than the 
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orbital period. 

If we assume that the entire trajectory of a spot is contained on the part of the primary 
star eclipsed by the secondary star then we can estimate the obliquity of the system (18) to 
be smaller than arctan(i?s/i?A) ~ 20°. The obliquity is likely to be smaller since we have 
detected more than 10 spots receding with different velocities, and these different velocities 
could be due to spots at different latitudes exhibiting differential rotation. We note that the 
obliquity of this target will be very hard to measure with the Rossiter-McLaughlin effect (40) 
due to its faintness, so additional investigation of its spots might be the preferred method to 
further constrain the obliquity. 

In principle, we can use the spectroscopic V TOt sin i together with an estimated rotation pe- 
riod and size of the primary star to obtain information on the inclination of the primary star. The 
spectroscopic observed V mt smi = 4.1 ± 0.5 km s -1 , while the inferred V rot = 2nR,A/P m t = 
6.3 ± 0.2 km s -1 . This would imply a highly inclined star (i s ps 40°). Note, however, that 
the measured value of the rotational velocity is below the resolution of the spectra, so its value 
should be treated with caution. In addition, differential rotation can make it harder to compare 
the surface integrated projected rotational velocity V rot sin i to the equatorial rotational velocity 
2nR A /P rot (39). 

1.8 Transit times for the inner and outer planet and the search for addi- 
tional transits 

All of the transits of the outer planet and about half of the transits of the inner planet are evi- 
dent in the SAP light curves before any detrending. The rest became visible when the data are 
carefully detrended. A symmetric polynomial "U-function" template with an adjustable width 
and depth was used to estimate the times of mid-transit and their durations. A cubic polynomial 
was used to detrend each transit using five different duration windows around the transit, and the 
best-fitting one was adopted. The fits for each transit were iterated to determine the best-fitting 
time of mid-transit and the duration. This method worked well in some cases and failed to 
converge in other cases. In cases where the convergence failed, the time was estimated using an 
interactive plotting program, and an uncertainty of 30, 60 or 90 minutes was assigned based on 
the judged quality of the transit. Table S4 gives the measured times and durations and their un- 
certainties, and the corresponding model times and durations. We note that the measured times 
and the durations were only used to establish starting models for the photometric-dynamical 
models described below. The actual (detrended) light curve was modeled directly. 

One "orphan" transit occurring about 12 hours after a transit of planet b was noticed in 
the Q12 data (Fig. IS 1 51 ). This transit cannot be accounted for by either the inner or the outer 
planet (the intervals between the nearest transits are 0.5 days and 127 days, respectively). To 
estimate the significance of the orphan event, a model consisting of two Gaussians was fit to the 
segment of the detrended light curve shown in Fig. IS 1 51 which contains 103 data points. The 
uncertainties on each point were scaled to give xt = 1 f° r 96 degrees of freedom. The Gaussian 
in the model at the location of the orphan was replaced by the background level of 1 .0 and the 
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resulting x 2 value increased to 205.5, giving a formal significance of ~ 10.5a. 

No other orphan transits with a significance of > 3a were found using visual searches. An 
automated search algorithm, dubbed the "Quasi-periodic Automated Transit Search" (QATS) 
was also used to search for additional transits. The QATS algorithm can allow for unequal time 
intervals between the transit events. For a given trial period for a potential planet, the expected 
transit duration at each time in the light curve is computed using a circular orbit for the planet. 
The data are corrected for the different transit durations and shifted to a common phase to 
increase the signal-to-noise (the correction for the different widths is quite good, provided the 
planet's orbit is nearly circular). A "periodogram" is constructed by plotting the significance 
versus the trial period. QATS detected the inner planet at high significance. Unfortunately 
QATS is very sensitive to detrending errors for longer periods, and in fact did not detect the 
outer planet. No additional planets with periods shorter than 150 days were detected at the 
significance level of the inner planet. 

Although the overall duty cycle of the data collection by Kepler is quite high, a non-trivial 
amount of the light curve is occupied by the stellar eclipses, which in the case of Kepler-47 is 
3 — 6 times more than it is for Kepler-16, 34, 35, and 38. Although one could in principle find 
transits during primary and secondary eclipses, in Kepler-47 this is extremely difficult owing 
to the effects of star-spots. The primary and secondary eclipse durations together are 0.014 in 
orbital phase, which is 0.104 days. A combined total of 256 primary and secondary eclipses 
were observed, giving a total of 26.62 days lost for the purposes of transit searches, lowering 
the duty cycle to 83.8%. 

Finally, if the transit is due to another planet in the Kepler-47 system, its radius would be 

4.5 Earth radii. Without more transit events, it is nearly impossible to determine what the 
orbital period of such a planet would be. If its orbit is more inclined relative to the other planets, 
it would not necessarily transit the stars near each conjunction. In addition, if there is precession 
of the orbit, it is possible for sequences of transits to come and go over long time scales. Thus, 
the orphan transit could in principle belong to a planet in between the inner and outer one, in 
spite of the lack of other observed transits. 

1.9 Photometric-dynamical model 

We modeled the Kepler light curve of Kepler-47 using a dynamical model to predict the motions 
of the planets and stars, and a eclipse/transit model to predict the light curve. 

1.9.1 Description of the model 

The "photometric-dynamical model" refers to the model that was used to fit the Kepler photom- 
etry. This model is analogous to that described in the analyses of KOI- 126 (13), Kepler-16 (6), 
Kepler-34 and Kepler-35 (7), Kepler-36 (41), and Kepler-38 (8). 

Four bodies were involved in this problem; however, the planets' gravitational interaction 
with the stars and with each other was determined to be observationally negligible. We therefore 
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assumed the planets to be massless in our model. The motion of the stellar binary was Keple- 
rian and could be predicted analytically. The planets were modeled as orbiting in the two-body 
potential of the stars. The motion of each planet was determined via a three-body numerical 
integration. This integration utilized a hierarchical (or Jacobian) coordinate system. In this 
system, r b (r c ) is the position of Planet b (Planet c) relative to the center-of-mass of the stellar 
binary (which corresponds to the barycenter in this approximation), and r EB is the position of 
Star B relative to Star A. The computations are performed in a Cartesian system, although it 
is convenient to express r b (r c ) and r EB and their time derivatives in terms of osculating Kep- 
lerian orbital elements: instantaneous period, eccentricity, argument of pericenter, inclination, 
longitude of the ascending node, and time of barycentric conjunction: Pb, c ,EB, eb,c,EB, %,c,eb, 
^b,c,EB, &b,c,EB, T bjC:EB , respectively. We note that these parameters do not necessarily reflect 
observables in the light curve; the unique three-body effects make these parameters functions 
of time (and we refer to these coordinates as "osculating"). 

The accelerations of the three bodies are determined from Newton's equations of motion, 
which depend on r b (r c ), r EB and the masses (42, 43). For the purpose of reporting the masses 
and radii in Solar units, we assumed GM Snn = 2.959122 x 1(T 4 AU 3 day -2 and R Sun = 
0.00465116 AU. We used a Bulirsch-Stoer algorithm (44) to integrate the coupled first-order 
differential equations for r& ; EB and i^eb- 

The spatial coordinates of all four bodies at each observed time are calculated and used as 
inputs to model the light curve. The computed flux was the sum of the fluxes assigned to Star 
A, Star B, and a seasonal (being the four "seasons" of the Kepler spacecraft orientation) source 
of "third light," minus any missing flux due to eclipses or transits (only planetary transits across 
Star A were computed, those across Star B are not significant in the Kepler data). The loss 
of light due to eclipses was calculated as follows. All objects were assumed to be spherical. 
The sum of the fluxes of Star A and Star B was normalized to unity and the flux of Star B was 
specified relative to that of Star A. The radial brightness profiles of Star A and Star B were 
modeled with a linear limb-darkening law, i.e., I(r)/I(0) = 1 — ui(l — \/l — r 2 ) where r is the 
projected distance from the center of a given star, normalized to its radius, and u is the linear 
limb-darkening parameter. The limb darkening coefficient of Star B was fixed (to u = 0.5); 
letting it vary freely resulted in a negligible change to final parameter posterior. 

The radial velocity of Star A was computed from the time derivative of the position of Star 
A along the line of sight (analytically, in this case) and compared to the radial velocity data. 

The continuous model is integrated over a 29.4 minutes interval centered on each long ca- 
dence sample before being compared to the long cadence Kepler data. 

1.9.2 Local detrending of Kepler data 

The Kepler light curve ("SAP_FLUX" from the standard fits product) for Kepler-47, spanning 
twelve Quarters, is reduced to only those data within 0.5-1 day of any primary or secondary 
eclipse or any transit of either planet. As noted above, some data are missing as a result of 
observation breaks during Quarterly data transfers or spacecraft safe modes. 
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Each continuous segment of data has a local cubic correction in time divided into it. The 
parameters of this polynomial correction are found through an iterative process, as described as 
follows. In the first step, we masked the eclipses of the stars and the transits of the planets and 
then performed a robust nonlinear least-squares fit to each continuous segment. The data, having 
divided out this correction, were then "fit" with the photo-dynamical model by determining the 
highest likelihood solution from a Markov Chain Monte Carlo simulation. The best-fit model 
was then divided into the data and the local nonlinear fits were recomputed (this time without 
masking the eclipses and transits). This process was repeated until the corrections converged to 
a sufficient tolerance. 



1.9.3 Specification of parameters 

A reference epoch for the three-body integration was specified for each planet near a particu- 
lar transit. Those epochs were chosen to be 2,454,969.216 BJD and 2,455,246.6545 BJD for 
planets b and c, respectively. 

The model has 33 adjustable parameters. Two parameters are related to Star A: the stellar 
density times the gravitational constant, GpA, and the stellar mass times G, GMa- One param- 
eter gives the mass ratio of the stars, q = M B /M A . Six parameters encode the eccentricities 
and arguments of pericenter for the planetary and stellar orbits about the barycenter in a way 
that reduces nonlinear correlations: 

h, c = y/e^sinu btC (1) 
h,c = y/eb^cosu byC (2) 
H = e EB smu E B (3) 
K = e EB cos uj eb (4) 

The remaining osculating parameters, 11 in total, are the periods P bc , Peb, the orbital 
inclinations i b;C , Ieb, the times of conjunction with barycenter, T b>c , Teb and the difference 
between the nodal longitudes of the planets to that of the stellar binary AQ b>c . The absolute 
nodal angle relative to North of the stellar binary cannot be determined and was fixed to zero in 
practice. 

Three more parameters are the relative radii of Star B and the planets to that of Star A: 
r B = R b /Ra and r b>c = R b , c / Ra- One parameter, u, parameterizes the linear limb darkening 
law for Star A (described above). Another parameter gives the relative flux contribution of Star 

B, F b /Fa- 

A single parameter, a LC , describes the width of the probability distribution for the photo- 
metric noise of the long cadence observations, assumed to be stationary, white and Gaussian- 
distributed. 

Three parameters characterize the radial velocity measurements: the constant offset of the 
radial velocity, 7, the offset between the HET and HJST velocities, A7, and a "stellar jitter" 
term, a R v, which contributes to the measured errors for each radial velocity observation, in 
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quadrature. Because only one Keck observation was made, this radial velocity could not be 
offset to match the HET and HJST velocities in a sensible way, and therefore was omitted in the 
modeling. 

Additionally, we specify 4 more parameters describing the relative extra flux summed in 
the aperture. The four parameters specify the constant extra flux in each Kepler "season." The 
Kepler spacecraft is in one of four orientations during a year; a constant level of "third light" is 
assumed for all Quarters sharing a common season. 

1.9.4 Priors and likelihood 

We assumed uniform priors for all 33 parameters. For the eccentricity parameters, this corre- 
sponds to uniform priors in eb, c and uib,c, but a prior that scales as ceb for the stellar eccentricity. 
This eccentricity is sufficiently determined that this non-uniform prior does not dominate the 
posterior distribution. 

The likelihood £ of a given set of parameters was taken to be the product of likelihoods 
based on the photometric data and radial velocity data: 



where AF^ is the ith photometric data residual, olc is the width parameter describing the 
photometric noise of the long cadence data, AVj is the jth radial velocity residual, Oj is the 
uncertainty in the jth radial velocity measurement and a^v is the stellar jitter term added in 
quadrature with the <jj. 

1.9.5 Best-fit model 

We determined the best-fit model by maximizing the likelihood. The maximum likelihood solu- 
tion was found by finding the highest likelihood in a large draw from the posterior as simulated 
with a Markov Chain Monte Carlo simulation as described below. Fig. IS 161 shows 18 transits 
of the inner planet, and Fig. IS 171 shows the residuals (observed data minus the model). With a 
few exceptions, there are no strong patterns in the residuals. Fig. IS 181 shows the model fits and 
the residuals for the outer planet. There are no patterns evident in the residuals. The residuals 
for the fits to the primary eclipses are shown in Fig. IS 191 and the residuals for the fits to the 
secondary eclipses are shown in Fig. IS201 Spot crossing events are evident in many of the pri- 
mary eclipses. Fig. IS2 II shows the radial velocity measurements and the best-fitting model and 
the residuals of the fit. Generally the absolute value of the radial velocity residuals is less than 
about 200 m s _1 . 



£ oc a L c LC exp - ^ 




(5) 
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The photometric noise parameter, cr^c, has a best-fit value of glc = 629.5 ppm. For com- 
parison, the root-mean- square deviation of the best-fit residuals is 626.9 ppm. This is similar 
to the expected noise in the light curve of 635 ppm as estimated using an on-line tool provided 
by the Kepler Guest Observer OfficeQ, where we used an apparent magnitude of Kp=15.18 and 
20 pixels in the aperture. For this ctlc, the x 2 -metric for the photometric data is x 2 — 10576 
with 10629 degrees of freedom. If we fail to include planet b in our model (by setting its ra- 
dius to zero), the x 2 increases by A% 2 = 343.4. If we ignore planet c, the x 2 increases by 
A X 2 = 248.2. 

The stellar jitter parameter, <Jrv, has a best-fit value of grv = 0.31 km s _1 . The value of 
X 2 for the radial velocity data alone is x 2 — 8.85 for the 10 radial velocity observations. 

Fig. IS22I shows schematic diagrams of the Kepler-47 orbits. The projected orbits of planets 
b and c cross the projected disk of the primary, and so transits of both planets across the primary 
occur, as do occultations of both planets by star A. The former events are observed, whereas 
the latter events are not observable given the noise level. On the other hand, owing to its 
small radius, the projected disk of star B does not intersect the projected orbits of the planets, 
and as such no transits of star B or occultations due to star B occur for the best-fitting orbital 
configuration. Due to the uncertainties in the relative nodal angles, transits of the planets across 
star B might occur for a subset of the acceptable solutions. However, even if transits across star 
B did occur, the expected transit depth would be ~ 30 times weaker than the transits across the 
primary, and would not be observable in the light curve given the noise level. 

1.9.6 Parameter estimation methodology 

We explored the parameter space and estimated the posterior parameter distribution with a Dif- 
ferential Evolution Markov Chain Monte Carlo (DE-MCMC) algorithm (45). 

We generated a population of 100 chains and evolved them through approximately 200,000 
generations. The initial parameter states of the 100 chains were randomly selected from an over- 
dispersed region in parameter space bounding the final posterior distribution. The first 10% of 
the links in each individual Markov chain were clipped, and the resulting chains were concate- 
nated to form a single Markov chain, after having confirmed that each chain had converged 
according to the standard criteria. 

The parameter values and derived values reported in Tables S5 and S6 beside the best-fit 
values (see above), were found by computing the 15.8%, 50%, 84.2% levels of the cumulative 
distribution of the marginalized posterior for each parameter. Figure IS231 shows two-parameter 
joint distributions between all parameters. This figure is meant to highlight the qualitative 
features of the posterior as opposed to providing quantitative ranges. The numbers in that figure 
correspond to the model parameters in Table [S6] with the same number listed as in the first 
column, if available. 

Figures IS24I and IS25I show the posterior distribution in the eccentricity and argument of 
pericenter planes The distribution of the three-dimensional inclination between the planets' 

1 http://keplergo.arc. nasa.gov/CalibrationSN.shtml 
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orbits and the invariable plane is shown in Figure IS261 



1.9.7 Predicted ephemerides and transit parameters 

Tables S7 and S8 provide the predicted times of transit, impact parameters, normalized transit 
velocities and durations over 7 years, starting with Kepler Quarter Q13. 

1.10 ELC light curve models 

Although the secondary star is not detected spectroscopically, its temperature can be estimated 
using the temperature of the primary derived using SPC, and the temperature ratio derived 
from modeling the eclipses. In order to find the temperature ratio, and to have an independent 
check on the results from the photometric-dynamical model, we modeled the light and velocity 
curves using the Eclipsing Light Curve (ELC) code (46) with its genetic algorithm and Monte 
Carlo Markov Chain optimizers. The free parameters include the temperature ratio T b /Ta, the 
primary's limb darkening parameters x& and y& for the quadratic limb darkening law \I(p) = 
7 (1 — x(l — fi) — y(l — /i) 2 )], the orbital parameters (e, to, i), and the fractional radii Ra/cl 
and Rb/cl- The stellar masses and the orbital period were held fixed at the values derived from 
the photometric-dynamical model discussed above. 

In ELC, the shapes of the stars are computed using a "Roche" potential modified to account 
for nonsynchronous rotation and eccentric orbits (47,48). Given the mass ratio and the fractional 
radii, the volumes of each star are found by numerical integration. The effective radius of each 
star is taken to be the radius of a sphere with the same volume as the equipotential surface. In 
the case of Kepler-47, the stars are very nearly spherical. For the primary at periastron, the ratio 
of the polar radius to its effective radius is 0.99988, and the ratio of the radius along the line 
of centers to the effective radius is 1.00007. The amplitude of the out-of-eclipse modulation in 
the light curve due to ellipsoidal variations, reflection, and Doppler boosting is on the order of 
400 ppm, which is ps 75 times smaller than the modulation due to star-spots. Thus the use of 
spherical stars in the photometric-dynamical model is a very good approximation. 

Since the numerical integrations are very CPU intensive, ELC has a fast "analytic" mode 
where the equations given in (49) are used. The normalized light curve was divided into 41 seg- 
ments containing two or three pairs of primary and secondary eclipses. These segments were 
modeled separately in order to help assess the systematic errors associated with the changing 
star-spots and the changes in the contamination from Quarter to Quarter. For each fitting pa- 
rameter, we computed the mean of the best-fitting values and the standard deviation. Table IS 1 01 
gives the mean values and standard deviations, which we adopt as la errors. 

Based on the temperature of the primary derived from the SPC analysis, and the temperature 
ratio found from the ELC models, we derive a temperature of 3357 ± 100 K for the secondary. 
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1.11 Upper limits on planetary masses 



Upper limits on the masses of the planets can be placed separately as follows. The mass of 
the inner planet is best constrained by the lack of eclipse timing variations due to gravitational 
perturbations from that planet. The planet will induce short-term eclipse timing variations with 
a period equal to the planet's period. It will also cause precession of the binary. Over the time- 
scale of a few years, the binary precession will cause a slight change in the phase difference 
between the primary and secondary eclipses, which can be observed as a slight difference be- 
tween the primary and secondary eclipse periods. Numerical simulations showed that in this 
case the stronger upper limit comes from the lack of short-term eclipse timing variations. A grid 
of masses for the inner planet was used, and equations of motion for a three body system were 
integrated, holding the orbital parameters of the binary at their best-fitting values (the nature 
of the perturbations on the binary are insensitive to anything except the planet's mass). The 
period and epoch of the binary was found, and the predicted times of eclipse were compared 
to the measured times. The x 2 value changes smoothly with the planet's mass, and going to 
X 2 = Xmin + 9 gives a 3a upper limit of 2.7 Jupiter masses for the uncorrected eclipse times 
and 2.0 Jupiter masses for the eclipse times corrected for the effect of star-spots. We adopt the 
latter value as the 3a upper limit on the mass of the inner planet. 

The mass of the outer planet was best constrained by light travel time (LTT) effects. We fit 
an LTT orbit to the corrected eclipse times, using a period of 303.13 days and constraining the 
eccentricity to be e < 0.2. While no convincing signal is seen at that period, the best- fitting 
orbit formally has a semiamplitude of 3.84 ± 1.84 seconds. Given the total mass of the binary 
and the period of the outer planet, we find a 3a upper mass limit of 28 Jupiter masses. 

1.12 Stability of orbits and limits on eccentricity 

We carried out an extensive study of the dynamics of the system and its long-term stability. 
The orbits of the two planets were integrated, numerically, for different values of their masses 
and orbital eccentricities. To determine an upper limit for the eccentricity of planet c, we held 
constant all orbital elements at their best-fit values and integrated the system varying the eccen- 
tricity of this planet. Results indicated that the system maintained stability for at least 100 Myr 
and for e c < 0.6. An examination of the semimajor axis, eccentricity, and orbital inclination of 
each planet during the course of the integrations showed that the variations of these quantities 
were negligibly small, supporting the idea that the two planets do not disturb each other's orbits. 
The results stayed unchanged when the masses of the two planets were increased to 0.21 and 
0.63 Jupiter masses, roughly ten times their plausible values based on the empirical mass-radius 
relations (14, 15). 

Both the photometric-dynamical model and stability simulations used a Newtonian 4-body 
numerical integrator. A more physical model would include the precession of the binary due to 
general relativity (GR) and the tidal and rotational bulges. Expressions for the rate of precession 
due to these effects (50) show that GR dominates, and it would cause a full periastron rotation 
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in ~ 6700 years. In the current observations, such precession would cause the period of primary 
and secondary eclipse to differ fractionally by ~ 10~ 7 , whereas the uncertainty of this quantity 
is 4.6 x 10~ 7 . The GR precession period is much longer than the periastron period of the 
planets - e.g., numerical integrations of planet c showing a ~ 560 year precession cycle due 
to the effective quadrupolar gravitational potential of the binary - so it has little dynamical 
importance. Therefore GR and other precession effects are neither detectable nor significantly 
change our assessment of stability, so GR has little dynamical importance. 

1.13 Comparison with stellar evolution models 

The reasonably precise absolute dimensions determined for the stars in Kepler-47 (4-5% rel- 
ative errors for the masses, and 1.8% for the radii) offer an opportunity to compare the mea- 
surements with models of stellar evolution. This is of particular interest for the late M-dwarf 
secondary, given that low-mass stars have shown discrepancies with theory in the sense that 
they are generally larger and cooler than predicted. These anomalies are believed to be due to 
stellar activity (51, 52). 

In Fig. IS27I we compare the measurements for the primary star with a stellar evolution 
track from the Yonsei-Yale series (53, 54), interpolated to the exact mass we measure. The 
metallicity of this model is set by our spectroscopic determination of [Fe/H] = —0.25, where 
we assume the iron abundance tracks the metallicity measurement from SPC. The model is 
consistent with the observations to within less than 2a, and the small difference may be due 
either to slightly biased spectroscopic parameters (temperature and metallicity) or a slightly 
overestimated mass for the primary star. As a check, we produced a photometric estimate 
of the temperature using available photometry from the KIC and empirical color-temperature 
calibrations along with the reddening listed in the KIC. The result suggests a value closer to 
5900 K than 5600 K, although we consider this evidence to be somewhat circumstantial and 
highly dependent on reddening. We confirmed that the level of agreement between theory 
and observation is independent of the adopted model physics by comparing the primary star 
parameters with BaSTI stellar evolutionary tracks (55), which yielded similar results as the 
Yonsei-Yale models. 

In Fig. IS28I we compare the measurements for both components against models from the 
Dartmouth series (56), which incorporate physical ingredients (equation of state, non-grey 
boundary conditions) more appropriate for low-mass stars. We find the radius of the secondary 
of Kepler-47 to be consistent with these models, which would be an exception to the general 
trend mentioned above, although the mass error is large enough (~4%) that the conclusion is 
not as strong as in other cases. Its temperature, however, is lower than predicted for a star of this 
mass by about 200 K. This deviation is in the same direction as seen for other low-mass stars. 
Because the secondary is so faint, we have no information on its activity level. Age estimates 
for the system from this figure and the previous one are somewhat conflicting, and only allow 
us to say that Kepler-47 is very roughly of solar age. 
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1.14 Details of habitable zone 



To determine the insolation limits of the habitable zone for Kepler-47 c, we follow the relations 
given by (24) that include the stellar temperature as well as luminosity. The temperature term 
accounts for the different relative amount of infrared flux to total flux, which is important for 
atmospheric heating. We use the criteria of a runaway greenhouse effect for the inner boundary 
and the maximum greenhouse effect for a cloud-free carbon dioxide atmosphere for the outer 
boundary. This is more conservative than the "recent Venus" and "early Mars" criteria, but less 
conservative than the "water loss" and "first carbon dioxide condensation" criteria (24). The 
secondary star emits only 1.7% as much energy as the primary star (and only 0.58% in the 
Kepler bandpass), so its contribution is neglected. The resulting insolation limits are shown as 
the dotted lines in Figure 3 (main text). The relations given in (57), i.e. a cloud-free atmosphere 
yield nearly identical limits. 

The average insolation for Kepler-47 c for a circular orbit is 87% of the Sun-Earth insolation, 
and varies by ~ 9% peak-to-peak. For an eccentricity of 0.2 the mean insolation is 89% and 
varies from 59% to 144% of the Sun-Earth value; for an eccentricity of 0.4, the mean is 96% 
and varies from 43% to 261%. Even in this latter case, which is ruled out at the 95% confidence 
level by the photometric-dynamical model, the mean is less than the Sun-Earth value, and it is 
the mean insolation that is most relevant for habitability (58). Thus for all allowed eccentricities, 
Kepler-47 c lies in the habitable zone. 

Because the primary star dominates the system both in luminosity and mass (so the primary 
star remains near the barycenter), the variation in insolation is relatively small for a circular 
planetary orbit. This is seen in the upper left panel of Figure 3, where the variations are caused 
by the 7.4-day orbit of the primary star. For large eccentricities, the variation in insolation are 
dominated by the non-circular orbit of the planet. 

It must be stressed that the habitable zone is defined such that liquid water could persist for 
a biologically significant time period on an Earth-like planet (i.e. with a terrestrial CO2/H2O/N2 
atmosphere, plate tectonics, etc.), and the formulations of (17, 24, 57) explicitly assume such 
conditions. For Kepler-47 c these conditions are not met. Nevertheless, the main point is that 
Kepler-47 c receives approximately the same amount of energy from its stars that the Earth 
receives from the Sun. 

While it neglects most atmospheric physics, the equilibrium temperature T eq of the planet 
is still a useful characterization. Assuming that the entire surface of the planet radiates isother- 
mally (i.e. the stellar insolation is efficiently advected around the planet), and for a Bond albedo 
of A B =0J, appropriate for a Neptune-size planet and 1 Sun-Earth insolation (59), a value of 
T eq ~ 200 K is found for eccentricities from 0.0 to 0.3. For A#=0.34, corresponding to the 
albedos of Jupiter and Saturn, T eq ~ 243 K. For an Earth-like albedo of 0.29, which is appropri- 
ate for a habitable-zone planet, T eq ~ 247 K. The greenhouse effect will lead to temperatures 
at the 1-bar pressure level that are higher by several tens of degrees. 
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Time (BJD-2,455,000) 

Fig. S4: SAP and detrended light curves. Top: The SAP light curves of Kepler-47 are shown. 
The colors denote the season and hence the spacecraft orientation where black is for Ql, Q5, 
and Q9, red is for Q2, Q6, and Q10, green is for Q3, Q7, and Ql 1, and blue is for Q4, Q8, and 
Q12. Bottom: The normalized and detrended light curve with the spot modulation removed is 
shown. Fifteen primary eclipses and thirteen secondary eclipses were missed during monthly 
data downloads, spacecraft rolls between Quarters, spacecraft safe modes, and interruptions 
caused by solar flares. 
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Time (BJD-2, 455,000) 

Fig. S5: Light curves showing spot modulations. The SAP light curves of Kepler-47 from 
Ql (top), Q5 (middle), and Q9 (bottom) are shown. The target appeared on the same detector 
module during these Quarters. A modulation in the out-of-eclipse regions due to star spots is 
evident. 
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Fig. S6: The autocorrelation function of the cleaned light curve with the eclipses removed. 



100p 

90- 
80- 
70- 
60- 
a 50- 

ra 
_i 

40- 
30- 
20- 
10- 



slope: 7.775 + 0.022 [day] 



6 

Peak number 



10 



12 



Fig. S7: The measured lag versus the peak number in the autocorrelation function. The 

measured lag of the peaks in the autocorrelation function displayed in Fig. |S6]is shown. The 
dashed line is a linear fit to these points. The slope of 7.775 ± 0.022 days is taken to be the 
rotation period of the primary star. 
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Fig. S8: Representative broadening functions. Broadening functions (BFs) from the 
HET+HRS (left) and the HJST+Tull spectrograph (right) are shown. The solid lines are the 
best-fitting Gaussians. The smaller peak in the HJST BF is due to the sky background. In all 
cases, the BF peak due to the sky was resolved from the object BF peak. 
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Fig. S9: Mean primary and secondary eclipse profiles. The observed profiles for the primary 
eclipse (dots, left panel) and the secondary eclipse (dots, right panel) arrived at after an iterative 
process. The Piecewise Cubic Hermit Spline (PCHS) models are shown as the solid lines. The 
increased scatter in the middle of the primary eclipse is most likely due to the effects of star 
spots. 
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Fig. S10: Observed minus computed curves for the stellar eclipses. Top: The Observed 
minus Computed (O-C) residual times of the primary eclipses. Coherent deviations of nearly 
two minutes are seen, with a quasiperiod of ~ 178 days. Middle: The O-C times of the primary 
after correction for the effects of star spots. No periodicities or trends are evident. Bottom: The 
O-C times for the secondary star. Note the change in the vertical scale. The error bars are not 
shown for clarity. The scatter is much larger, and no periodicities or trends are seen. 
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Fig. SI 1 : Lomb-Scargle periodograms of O-C curves. Top: A Lomb-Scargle periodogram 
of the O-Cs of the primary eclipses, before any corrections for star spots have been applied. 
The peak power occurs at a period of 179.2 days (dashed line). The expected beat period of 
w 170 days is indicated by the dotted line. Bottom: A Lomb-Scargle periodogram of the 
primary eclipse O-Cs after a correction for the effects of star spots has been applied. There is 
no significant power at any period. 
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Fig. S12: The correlation of residual O-C time and local slope near primary eclipse. The 

dependence of the primary eclipse O-C times on the local SAP light curve slope is shown. A 
clear correlation is seen. The best-fitting line has a coefficient of correlation of r = —0.80. 
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Fig. S13: The effect of star-spots on the primary eclipses. Upper Left: The observed eclipse 
light curves (black dots) for five consecutive primary eclipses are shown. A model with no 
spots (red curves) does not fit the data well, whereas a model with a spot that is occulted by the 
secondary star fits much better (blue curves). Upper Right: The residuals for the same eclipses 
are shown. As time passes (top to bottom) the residual feature from the no-spot model moves 
from the right side of the eclipse to the left. Lower: A section of the light curve spanning 
eclipse cycles 60-64 is shown. Notice how the local slope in the immediate vicinity of the 
primary eclipse slowly changes from cycle to cycle. 
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Fig. S14: Eclipse depth variation and its effect on the secondary star radius estimate. Top: 
The individual depths for each primary eclipse calculated with a one-spot model are shown (dif- 
ferent color correspond to different observing seasons). The depth changes with time because 
the fraction of the star covered by spots changes with time. There is also a hint that the depths 
change with the observing season (each season the star falls into a different CCD, changing 
the level of contamination). Bottom: A histogram of the inferred radius of the secondary star 
(black line) for each eclipse is shown. This demonstrates how the secondary star radius from 
the photometric-dynamical model (thick red line) is slightly underestimated (as expected), but 
the difference is not significant compared to the error bars on the measured radius (the 15.4% 
and 84.6% confidence levels are shown with dotted red lines). 
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Fig. S15: A segment of the Q12 light curve showing a transit of the inner planet and an 
orphan transit. An "orphan" transit that cannot be accounted for by the inner or outer planets 
appears near the middle of this data segment, about 12 hours after a transit of the inner planet 
(left). The solid line is a simple model consisting of two Gaussians used to find the mid-transit 
time of the orphan and to evaluate the significance of the event. 
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Fig. S 16: All observed transits of the inner planet. The complete set of planet b transits with 
the best-fitting model is shown. The color coding is the same as in Fig. [S4l 
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Fig. S17: The residuals of the model fits of the inner planet transits. The residuals of the 
model fits of the transits due to the inner planet displayed in Fig. IS16I are shown. The color 
coding is the same as in Fig. | 
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Fig. S18: The model fits and residuals of the transits of the outer planet. The model fits to 
the transits of the outer planet are displayed in the top panels, and the residuals are shown in the 
lower panels. The color coding is the same as in Fig.[ 
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Fig. S19: The residuals of the fits to the primary eclipses. The residuals during each primary 
eclipse are displayed. As expected, numerous spot crossing events are seen. 
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Fig. S20: The residuals of the fits to the secondary eclipses. The residuals during each 
secondary eclipse are shown. As these eclipses are total there is much less structure seen in the 
residuals, compared to the residuals for the primary eclipses displayed in Fig. IS 19[ 
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Fig. S21: The observed and model radial velocity curve for the primary. Top: The radial 
velocity measurements shown as a function of orbital phase and the best-fitting model. Bot- 
tom: The residuals of the fit. Measurements from each telescope+instrument combination are 
denoted with different symbols and colors. 
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Fig. S22: Schematic diagrams of the Kepler-47 orbits. Top: A face-on view of the stellar 
and planetary orbits found from the best-fitting model of the Kepler-47 system. The center of 
mass of the system is marked with the cross. The stars and the planets would not be seen at 
this scale, and so their positions are marked with boxes. Bottom: The view of the system as 
seen from Earth on an expanded scale is shown. The lines denote the projected orbits of the 
various bodies. Both planets can transit the primary star (labeled A). Transits of the secondary 
star (labeled b) are narrowly missed for the best-fitting orbital configuration. 
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Fig. S23: Two-parameter joint posterior distributions of primary model parameters. The 

densities are plotted logarithmically in order to elucidate the nature of the parameter correla- 
tions. The indices listed along the diagonal indicate which parameter is associated with the 
corresponding row and column. The parameter name corresponding to a given index is indi- 
cated in Table IS6l in the first column. 
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Fig. S24: Posterior distributions in eccentricity. 
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Fig. S25: Posterior distributions in the eccentricity and argument of pericenter planes. 
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Fig. S26: Posterior distributions in the inclination of the planetary orbits relative to the 
invariable plane. 
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Fig. S27: Comparison of the absolute dimensions of the primary of Kepler-47 against 
stellar evolution theory. The thick solid line shows an evolutionary track from the Yonsei- Yale 
series (53,54) interpolated to the measured mass of that star and its measured metallicity. The la 
uncertainty in the location of the track due to the mass error is indicated with the darker shaded 
area. The wider light shaded area includes the additional la contribution from the uncertainty 
in [Fe/H]. Isochrones from 1 to 13 Gyr (left to right) are shown with dashed lines. 
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Fig. S28: Isochrones in the mass-radius and mass-temperature planes. Isochrones from 
the Dartmouth models (56) corresponding to ages from 1 to 13 Gyr, compared against the 
measured masses, radii, and temperatures of the stars in Kepler-47. The oldest isochrone is 
indicated with a solid line, and the metallicity has been set to the spectroscopically determined 
value of [Fe/H]= —0.25. The error bars for the measurements are represented with the shaded 
boxes. Top: The mass-radius diagram. The inset shows an enlargement around the secondary, 
which appears to agree with the models. Bottom: The mass-temperature diagram, showing the 
secondary to be cooler than predicted. 
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Table S2: Radial velocities for Kepler-47. 



Date 


UT Time 


BJD 




telescope 


YYYY-MM-DD 




(2,455,000+) 


km s 




2012-04-10 


13:25:48.68 


1028.05942 


11.442 ± 0.011 


Keck 


2012-04-23 


09:11:27.36 


1040.90325 


33.534 ± 0.091 


HET 


2012-05-01 


09:52:55.08 


1048.93237 


35.458 ±0.171 


HJST 


2012-05-02 


07:23:45.95 


1049.82882 


24.430 ± 0.440 


HJST 


2012-05-04 


08:34:10.40 


1051.88474 


-21.957 ±0.159 


HJST 


2012-05-05 


08:08:26.83 


1052.86692 


-26.719 ±0.178 


HJST 


2012-05-06 


08:08:55.42 


1053.86729 


-9.150 ±0.122 


HJST 


2012-05-18 


07:35:21.15 


1065.83749 


-1.843 ±0.060 


HET 


2012-05-20 


07:37:07.78 


1067.83880 


-25.681 ±0.030 


HET 


2012-06-05 


06:29:24.15 


1083.79236 


-5.223 ±0.080 


HET 


2012-06-26 


08:03:52.39 


1104.85862 


-26.743 ± 0.086 


HJST 



Table S3: Spectroscopic parameters from SPC. 



parameter 


value 


T cS (K) 


5636 ± 100 


log g (cgs dex) 


4.42 ±0.10 


[m/H] (dex) 


-0.25 ±0.08 


V TOt sinz (km s _1 ) 


4.1 ±0.5 



Table S4: Times of stellar eclipses. 



cycle # 


primary 


corrected 


uncertainty 


cycle # 


secondary 


uncertainty 




time^ 


time 1 


(min) 




time 1 


(min) 


0.0 








0.4873910 


-33.12216 


2.18 


1.0 


-29.30630 


-29.30631 


0.39 


1.4873910 


-25.67634 


2.58 


2.0 


-21.85791 


-21.85777 


0.37 


2.4873910 


-18.23125 


2.18 


3.0 


-14.40955 


-14.40931 


0.42 


3.4873910 


-10.78077 


2.18 


4.0 


-6.96153 


-6.96144 


0.43 


4.4873910 


-3.33057 


2.18 


5.0 








5.4873910 


4.11649 


2.18 


6.0 


7.93529 


7.93560 


0.34 


6.4873910 


11.56631 


2.48 


7.0 








7.4873910 


19.01194 


2.28 


8.0 


22.83203 


22.83220 


0.34 


8.4873910 


26.46447 


2.38 



2 BJD-2,455,000 
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Table S4: (continued) 



Q 

y .vj 


10 9R050 


10 9R0R0 


11 


Q 4R71Q10 

y .'-to I j i vj 


11 Q1 1Q6 

JJ.y 1 J -7VJ 


9 OR 

Z..VJO 


10 


37 79RRQ 

J I . I LOO y 


17 79909 

J / . / £.y\J y 


19 


10 4R71910 


41 16049 


9 9R 


1 1 

1 1 ,\J 


45 17791 

" J . 1 / / Z 1 


45 17755 

T 1 J . 1 / / J J 


11 

\J . J J 


1 1 4R71910 


48 80785 

^O . O \J I O J 


9 1 R 

Z.. 1 O 


1 9 

1 z . VJ 


59 6954Q 


59 69570 


15 

VJ. J J 


19 4R71Q10 

IZ.-'rO / J ;7 1 VJ 






1 1 

J- J . VJ 


60 07494 

VJVJ.VJ / tZ't- 


60 07441 

VJVJ.VJ / ttJ 


15 

VJ. J J 


1 1 4R71Q10 

i j.'to / J ;7 1 VJ 






14 


67 5996R 


67 59976 

VJ / . JZ.Z. / VJ 


44 

VJ.'-r'-r 


14 4R71Q10 

1 H.tO / J ;7 1 VJ 


71 15005 


9 1 R 


1 5 

1 J .V 


74 Q70Q1 

/ *T. j7 / Vjj? 1 


74 Q70Q0 


41 

VJ.'-r 1 


1 5 4R71Q10 

1 J.'TO / J!71VJ 


7R 601 17 

/ O .VJVJ 11/ 


9 1R 

Z.JO 


1 6 

1 VJ.VJ 


R9 41 Q51 


R9 41 Q4Q 
oz.t i y-ry 


1Q 


16 4R71Q10 

1 VJ. TO / JZ71VJ 


R6 04R9Q 

OVJ.VJH'OZ.Z? 


9 1 R 


17 

1 / .VJ 


RQ R67Q5 
oy .ovj / j 


RQ R6R0Q 

O;? . OVJOVJ;? 


1Q 

vj. j y 


17 4R71Q10 

1 / .'"TO / J !7 1 VJ 


g^j 4QQQ7 


9 1 R 


1 R 

1 o . V/ 


Q7 1 1 647 

/ / Jl U"t / 


Q7 1 1 640 


97 

VJ.Z / 


1 R 4R71Q10 

1 O .tO / J!71VJ 


100 Q4R96 


9 1R 

Z.JO 


1 Q 

1 ;? . vj 


104 764R9 


1 04 76476 

1 V./'-r . / Ut / VJ 


9Q 

VJ.Z;7 


1Q 4R71Q10 

ly.TO / J 1 VJ 


1 0R 1Q407 

1 VJO . J ^'-tvj / 


9 1 R 

Z. 1 O 


90 


119 91116 

1 1 Z.Z 1 J1U 


1 1 9 91 9R5 

1 1 Z.Z 1 Z,OJ 


17 


90 4R71Q10 


1 1 5 R4101 


9 1 R 

Z.IO 


91 

Z 1 .V 


1 1 Q 661 5R 
1 1 y .vjvj i j o 


1 1 Q 661 1R 

1 1 17 . VJVJ 1 JO 


9R 

VJ.ZO 


91 4R71Q10 

Z.l.'-rO / J ;7 1 VJ 


1 91 9Q097 

1 Z. J .Z.Z7VJZ. / 


9 9R 

Z. ZO 


99 

ZZ.VJ 


197 10Q71 


197 10Q51 

1 / . 1 VJ37 J 1 


94 

VJ.Z1" 


99 4R71Q10 

zz.to / j y i vj 


1 10 715QQ 

ijvj. / ojyy 


9 1 R 

Z.IO 


91 

Z J .VJ 


1 14 55R1 6 


1 14 55R05 


1R 

VJ. J o 


91 4R71Q10 

Z.J.H-0 / J ;7 1 VJ 


1 1R 1 R757 

IJO.IO / J / 


9 9R 

Z. ZO 


94 

Z'T.VJ 


1 49 00644 

1 tz. . VJVJVJ'T'T' 


149 0061Q 
i'tz.vjvjvj j y 


19 

VJ. JZ, 


94 4R71Q10 

ZT.TO / J!71VJ 


145 61599 


9 1R 

Z.JO 


95 


14Q 45471 


14Q 45475 


19 


95 4R71Q10 


1 51 0R9R7 


9 9R 
z.zo 


96 

ZVJ.VJ 








96 4R71Q10 

ZVJ.H'O / J v 1 VJ 


1 60 514RR 

1 VJVJ. J J TOO 


9 4R 

Z .^tO 


97 

Z / .VJ 


1 64 151 56 

i un - . J J i j vj 


1 64 151 5R 

1 VJ'-r . J J i jo 


10 

VJ. JVJ 


97 4R71Q10 

Z / .'-+0 / J 1 VJ 


1 67 QR07Q 

i vj / .z70vj / y 


9 OR 

Z.VJO 


9R 
zo .vj 


171 R0006 


171 R0016 

1 / 1 . OVJVJ 1 VJ 


14 

VJ. JT 


9R 4R71Q10 

ZO ."O / JZ71VJ 


175 411 10 

1 / J .*"+ J 1 Jy 


9 1 R 

Z.IO 


9Q 

Z;7 .VJ 


17Q 94R07 


17Q 94R06 

1 / y .Z'-rOVJVJ 


95 

VJ.Z J 


9Q 4R71Q10 

Z.7.TO / J 1 VJ 






10 


1 R6 6Q696 

i ovj.vj;?vjzvj 


1 R6 6Q604 

i o vj . vj y vjvj'-r 


11 

VJ. J J 


10 4R71Q10 

JVJ. TO / J!71VJ 


1 Q0 19696 


9 1 R 

Z. 1 o 


11 

J 1 .V 


1Q4 14450 

1 y~r. 1 t-'t- J VJ 


1Q4 144R6 

i y-x . i T-to vj 


11 

VJ. J J 


11 4R71Q10 

Jl.TO / J!71VJ 


1 Q7 77591 

ly / ./ / JZ. J 


9 9R 

Z.ZO 


19 


901 5Q977 

z.v_/ i . j z / I 


901 5Q107 

zi_/ i . j y jvj I 


40 

VJ.T-VJ 


19 4R71Q10 

JZ,.TO / J!71VJ 


905 99579 

ZUJ . Z. Z. J / z. 


9 1 R 

Z.IO 


11 

J J .VJ 


90Q 041 9Q 

ZVJ;7 .VJ^T 1 Z.7 


90Q 041 61 

zjyjy .vj'-r i vj i 


4R 

VJ.'-rO 


11 4R71Q10 

JJ-TO / J!71VJ 


919 67179 


9 1 R 

Z.IO 


14 








14 4R71Q10 

J'r.TO / J!71VJ 


990 1 1Q11 

Z.Z.VJ. 1 17JJ 


9 1 R 

Z.IO 


15 

J J .V 


991 Q1R96 

Z^Z^J .7JOZU 


991 Q1R41 


41 

VJ.'-T 1 


15 4R71Q10 

J J . TO / J!71VJ 


997 567R1 

Z.Z. / .JVJ / O J 


9 9R 

Z.ZO 


16 

J VJ.VJ 








16 4R71Q10 

JVJ. tO / JZ71VJ 


915 0171R 

Z/JJ .VJ 1 / J o 


9 1 R 

Z. 1 o 


17 

J / .VJ 


91R R1551 


91R R1546 

Z-JO.OJ Jt-U 


11 

VJ. J 1 


17 4R71Q10 

J / .'"TO / JZ?1VJ 


949 46451 

Z-T-Z. . t TVJ t T J 1 


9 1 R 

Z.IO 


38.0 


246.28389 


246.28401 


0.33 


38.4873910 


249.91181 


2.18 


39.0 


253.73230 


253.73247 


0.29 


39.4873910 


257.36186 


2.58 


40.0 


261.18019 


261.18031 


0.25 


40.4873910 


264.80933 


2.18 


41.0 


268.62872 


268.62889 


0.27 


41.4873910 


272.25848 


2.18 


42.0 








42.4873910 


279.70576 


2.18 


43.0 


283.52551 


283.52545 


0.33 


43.4873910 


287.15579 


2.58 


44.0 


290.97418 


290.97397 


0.32 


44.4873910 


294.60501 


2.18 


45.0 


298.42259 


298.42234 


0.38 


45.4873910 


302.05267 


2.18 


46.0 


305.87113 


305.87077 


0.28 


46.4873910 






47.0 


313.31965 


313.31911 


0.33 


47.4873910 


316.94980 


2.18 
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Table S4: (continued) 



48 


190 76801 

/ VJOVJJ 


120 76741 

JZ.VJ. / VJ / t 1 


11 

VJ. J 1 


48 4871910 

l tO.'tO / Jy 1 VJ 


124 19888 

JZ.t. Jy OOO 


2 1 8 

Z.. 1 O 


49 

t y . \J 


128 21615 


128 21 580 


11 

\J . J J 


49 4871910 


111 84676 

JJ1 .OtVJ / VJ 


2 08 


50 

J vy . vy 


115 66454 

J J J .UUtJt 


115 66412 

J J J . VJ VJ T J Z. 


10 

VJ . J VJ 


50 4871910 

J vy . T O / Jy 1U 


HQ 29109 


2 28 


51 


141 1 1 269 


141 1 1 255 


27 

VJ . z. / 


51 4871910 

JI.tO / Jy IV/ 


146 74012 

JtVJ. / tVJJZ. 


2 1 8 


52 


150 56108 

J JU .JVJ 1 VJO 


150 56100 

J J VJ . J VJ 1 VJVJ 


17 

VJ. J / 


52 4871910 

JZ,.tO / J7 lU 


154 1 8991 

J Jt . 1 Oy y 1 


2 28 


51 

J J .VJ 


158 0091 6 

JJO .UU7 1 VJ 


158 00920 

J J O . \J\J y Z. V / 


17 

VJ. J / 


51 4871910 

JJ-TO / J7 IV/ 


161 61907 
jvj 1 ,\joy\j 1 


2 1 8 


54 


165 45750 


165 45727 

J VJ J ,*T J / Z. / 


15 

VJ. J J 


54 4871910 

Jt .to / Jy lu 


169 08855 

JVJy.VJOOJJ 


2 1 8 


J J .VJ 


172 90581 

J / z..yvj j o J 


172 90551 


17 

VJ. J / 


55 4871910 

JJ.tO / J7 lU 


176 51671 

J / VJ. J JVJ / J 


2 18 


56 
j vj. vj 


180 15414 


180 15414 


14 

VJ. JT- 


56 4871910 

JU.tO / JylVJ 


181 98451 

J O J . y Ot J J 


2 18 


57 

J / .VJ 


187 80274 

JO / .OVJZ. / i 


187 80104 

JO 1 ■ OVJ JVJt 


29 

VJ.Z.y 


57 4871910 

J / .'-+0 / Jyivj 


191 41101 

J y 1 .tj Jul 


2 18 


58 

J O . VJ 


195 25105 

u 1 VJ J 


195 251 58 

J J? J .Z. J 1 JO 


28 

VJ.Z.O 


58 4871910 

JO.tO / JylVJ 


198 881 91 

J y O .OO 1 y J 


2 28 


5Q 


402 69912 
L r\j z. .yjy y j 


402 69981 


19 


59 4871910 


406 11048 


2 28 
_. ._.o 


60 

VJVJ.VJ 


410 14751 

1 1 VJ. 1 t / J J 


410 14797 

t IVJ. 1 T- / 7 / 


19 

vj. j y 


60 4871910 

VJVJ. to / Jyivj 


41 1 77876 

tlj./ / O/ VJ 


2 48 

Z..TO 


61 


417 59598 


417 59649 

t 1 / . J y VJty 


17 

VJ. J / 


61 4871910 

VJ 1 . to / Jyivj 


421 22705 

TZ. 1 / VJ J 


2 28 


VJZ..VJ 


425 04448 

tZ. J .VJtttO 


425 04486 

tZ. J .VJttOVJ 


41 

VJ.t 1 


62 4871910 

UZ,.tO / JylVJ 


428 67664 

tZ-O.VJ / VJVJt 


2 1 8 


61 

VJ J . 1/ 


412 49299 


412 49115 


19 

vj. j y 


61 4871910 

UJ-tO / JylVJ 


416 12677 

t J VJ. 1 Z.VJ / / 


2 48 

Z..TO 




419 941 60 


419 941 59 


18 


64 4871910 


441 57189 

L r L rJ .J / Joy 


2 28 


65 

UJ .V/ 


447 18994 


447 18962 


41 


65 4871910 
uj.ho / jy ivj 


451 01777 


2 28 


66 

VJVJ.VJ 


454 81858 

L rJ L r.OJOJO 


454 81824 


42 

VJ.tZ. 


66 4871910 

VJVJ . t O / ^jy 1\J 


458 46971 

t J O .tVJy / J 


2 28 


fn 








67 4871910 

VJ / .tO / JylVJ 


465 91 872 

tVJ J .y 1 O / Z. 


2 1 8 


68 

VJO . VJ 


469 71590 

tvjy . / *jjy\j 


469 7151 1 

tVJy • / J J 1 1 


49 


68 4871910 

VJO. to / JylVJ 


471 16671 

t / J . JVJVJ / J 


2 1 8 


69 

vjy .vj 


Ail i 8426 

t / / . 1 OtZ.VJ 


477 1 8114 

t / / . 1 OJJ L r 


41 


69 4871910 

VJy. to / JylVJ 


480 80981 

tOVJ.OVJyOJ 


2 1 8 

z.. 1 


70 

/ VJ.VJ 


484 61271 

tOt.VJ JZ. / J 


484 611 84 

tOt. VJ J lOt 


42 

VJ.tZ. 


70 4871910 

/vj.to / Jyivj 


488 26171 

tOO.ZUJ / 1 


2 1 8 

Z. . 1 


71 


492 08010 


409 07982 

t y ZZ . VJ / yOZ. 


29 

VJ.Z.y 


71 4871910 

/ 1 .to / Jy iv 


495 71 271 

tyJ./ 1Z / J 


2 08 

Z..VJO 


72 

/ Z..VJ 


499 52841 


499 52811 


15 

VJ. J J 


72 4871910 

/ Z..tO / Jyivj 


501 1 5841 

JVJ J . 1 J ot 1 


2 28 

Z. . Z. O 


71 

/ J .VJ 


506 97666 

JV/VJ.!? / VJVJVJ 


506 97688 
u\j\j.y 1 vjoo 


17 

VJ. J / 


71 4871910 

/ J. to / Jyivj 


510 60646 

J IVJ.VJVJVJtVJ 


2 28 


74 

/ t .VJ 


514 42490 


514 42526 


11 

VJ. J J 


74 4871910 

/t.to / Jyivj 


518 05528 

J 1 O .VJ J JZ.O 


2 1 8 

Z. . 1 O 


75 


521 87114 


521 87156 


11 

VJ. J 1 


75 4871910 


525 50508 


2 28 


If, 

1 VJ.VJ 


529 121 82 


529 1221 6 

j z.y . j z.z* 1 vj 


28 

VJ.Z.O 


76 4871910 

/VJ.tO / Jyivj 


512 95444 

J JZ..y Jttt 


2 18 

Z.. J O 


77.0 


536.77016 


536.77052 


0.27 


77.4873910 


540.40058 


2.28 


78.0 








78.4873910 


547.84725 


2.38 


79.0 


551.66689 


551.66678 


0.32 


79.4873910 






80.0 








80.4873910 






81.0 








81.4873910 


570.19645 


2.88 


82.0 


574.01232 


574.01196 


0.31 


82.4873910 


577.64083 


2.38 


83.0 


581.46071 


581.46061 


0.32 


83.4873910 


585.09178 


2.38 


84.0 


588.90872 


588.90902 


0.35 


84.4873910 


592.54015 


2.18 


85.0 








85.4873910 


599.98822 


2.08 


86.0 


603.80549 


603.80542 


0.56 


86.4873910 


607.43557 


2.28 
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Table S4: (continued) 



R7 


61 1 95410 


61 1 95400 

VJ 1 1 . Z..J H-VJVJ 


51 

VJ. u J 


R7 4R71910 

O / . H" O / J ;? 1 VJ 


614 RR579 

VJ It.OOJ / Z. 


9 OR 

Z..VJO 


RR 


61 R 70957 

VJ _I_ O . / \JZ^.J / 


61 R 70941 

VJ _I_ O . / Z.^ 1 


48 


RR 4R71910 


699 331 13 


9 1R 

Z. . J o 


RQ 


696 1 5074 


696 1 5090 

VJZ.VJ. 1 JUyU 


54 


R9 4R71910 


699 7R096 


9 1 R 

Z.. 1 o 


90 


fill 5QQ1Q 


fi^ 5QQ59 


47 

VJ.'-T / 


90 4R71910 






91 








91 4R71910 

yl.tO / ,J;7 1 VJ 


644 67R01 

Utt . VJ / OVJ J 


9 1 R 

Z. . J_ o 


99 


64R 49691 

VJ'-rO .tyUZ- J 


64R 4959R 


97 

VJ.Z. / 


Q9 4R71Q10 

yz.TO / oy 1 VJ 


659 19511 

VJJZ.. 1 Z. J J J 


9 1 R 

Z.. _I_ o 


91 


655 9445R 


655 94414 

VJ,J,J . ytt 1 i 


17 

VJ. J / 


91 4R71910 

y J-t-O / J ;? 1 VJ 


659 57495 

vjj y .j / t-z.-j 


9 1 R 

Z. . 1 o 


94 


661 1911 R 


661 19964 

VJVJ J . J ^'Z.VJ'-r 


11 

VJ. J 1 


Q4 4R71910 

/t.tO / J ;7 1 VJ 


667 09407 

VJVJ / .VJZ. L +VJ / 


1 9R 

l.yO 


95 


670 R4174 


670 R4141 

VJ / VJ. O^T 1 *T J 


19 


95 4R71910 

yJ-T-O / J ;? 1 VJ 


674 4701R 

VJ / T-.t / VJ J o 


9 1 R 

Z. . 1 O 


96 
-7vj. vj 








96 4R71910 

!7VJ. t TO / J!71VJ 


6R1 99000 
vjo i . y 


9 9R 

Z. . Z.O 


97 
/ .vj 


6R5 71R07 

VJO J . / JOVJ / 


6R5 71795 

vjo J . / J / y ,j 


19 

VJ. JZ 


97 4R71910 

y I .ho / j ;? i vj 


6R9 17001 

VJO-7.J / VJVJ J 


9 9R 

Z. . Z.O 


QR 

7U.U 


6Q1 1 R69Q 


6Q1 1 R690 


9Q 


QR 4R71Q10 


6Q6 R1 5R1 

U7U.O UOl 


9 OR 


QQ 
yy.\j 


700 61456 


700 6145R 


9R 

V7.Z.O 


QQ 4R71Q10 


704 96500 


9 OR 


100 


70R 0R104 

/ VJO .VJO JUT" 


70R 0R105 

/ VJO .VJO J VJ J 


99 

VJ.Z.;? 


100 4R71Q10 

ivjvj.'to / oy ivj 


71 1 71 161 

/ 11./ 1JU1 


9 OR 

Z..VJO 


101 

1 I J 1 . 1./ 


71 5 511 1 R 

/ 1J.JJ1 io 


71 5 511 10 

/ 1J.JJ 1 JW 


9R 

VJ.Z.O 


101 4R71Q10 

i vj i .t-o / oy ivj 


71 Q 16004 

/ iy . l vjvjvj't 


9 1 R 

Z.. _I_ o 


10? 

i vjz..vj 


799 Q7Q14 

/ Z.Z..;? I 7JT 


799 Q7Q55 


99 


109 4R71Q10 

lVJZ..'-rO / J;7±VJ 


796 60R75 

/ Z.VJ.VJVJO / J 


9 1 R 

Z.. J_ o 


101 


7^0 4.977R 


710 49R05 


99 

VJ.Z.Z. 


101 4R71Q10 


714 05R0Q 


9 1 R 

Z.IO 


104 

i vj'-r .vj 


7^7 R7697 

/ J / .O / VJZ. / 


717 R7690 

/ J / .O / VJZ.VJ 


97 

VJ.Z. / 


104 4R71Q10 

IVJ'-t.'-tO / J IVJ 


741 50754 


9 1 R 

Z.. 1 o 


105 


745 19471 

/ T^«J .JLT / .J 


745 19479 

/ '"T.J . J Z.'T" / Z. 


19 

VJ. jz 


105 4R71Q10 

i vjj .'-to / jy i\j 


74R Q55Q6 

/ to .yj j yyj 


9 1 R 

Z.. 1 o 


106 

1 VJVJ.VJ 


75? 7711 9 

/ jz.. / / j> Ly 


759 7711 1 

/ ,J Z. . / / J 1 1 


97 

VJ.Z. / 


106 4R71Q10 

1 VJVJ. tO / Jy IVJ 


756 40446 

/ JVJ.'~rVJ t -r'-rVJ 


9 9R 

Z.. Z.O 


107 

1 VJ 1 , VJ 


760 991 5R 

/ VJ I J . Z. Z. 1.JO 


760 991 69 

/ VJ V J . Z. Z. 1 uz 


10 

VJ. JVJ 


107 4R71Q10 

ivj / .to / oy i vj 


761 R51 1 6 

/ VJJ .OJ 1 1 VJ 


9 9R 

Z.. Z.O 


10R 


767 669R6 

/ VJ / .VJVJZ70VJ 


767 66997 

/ VJ / . VJVJZ/v / 


19 

VJ. JZ. 


10R 4R71Q10 

1 VJO .'"TO / JZ71VJ 


771 1006Q 

/ / 1 .JVJVJVJ^' 


9 9R 

Z.. Z.O 


109 


775 1 1 R1 R 


775 1 1 R94 

/ / J . 1 1 OZt 


11 

VJ. J J 


10Q 4R71Q10 

i xjy , t ~ro i oy ivj 


77R 74507 

/ / O . / *T J VJ / 


9 9R 

Z. . Z. O 


110 

1 1 VJ.VJ 


7R9 56649 


7R9 56650 

/ OZ.. JVJVJJVJ 


47 

VJ.'-T / 


1 10 4R71Q10 

i ivj.'-ro / oy ivj 


7R6 1QR17 

/ OVJ. 1 yOJ / 


9 4R 


1110 

J- 1 J- .V 


790 01 51 R 


790 01 591 


15 

VJ. J J 


1 1 1 4R71Q10 

_I_ 1 _I_ .1"0 / Jy J_ VJ 


7Q1 64R04 

/ y j .vj'-tovjh' 


9 OR 

Z..VJO 


119 


797 46174 

i y i . t rvj,j / t" 


797 461R7 

/ !7 / .'-rVJ J O / 


17 

VJ. J / 


1 19 4R71Q10 

J. IZ.tO / Jy J_ VJ 


R01 0Q090 

OVJ 1 . vj y VJ Z. 1 J 


9 9R 

. z. o 


1110 


R04 91 944 

OVJ'-T . -7 1 Z.^T'T - 


R04 91 911 

OVJ'-r. y 1 Z. J 1 


4R 

VJ.^+O 


1 1 1 4R71Q10 


R0R 54909 

OVJO . J'-pZ.VJZ. 


9 1 R 

Z.. _I_ o 


114 


«19 16061 

O 1 Z^. JVJVJVJJ 


R19 16010 

O 1 Z.. J VJVJJVJ 


99 

VJ.Z.;? 


1 14 4R71Q10 
i i t.ho / oy ivj 


R1 5 QQ09Q 

o u .yy\jz.y 


9 1R 


115 


R1 9 R0R79 

O 1 y.OUO / y 


R1 9 R0R61 

Oly. OVJOVJ 1 


90 

VJ.Z.VJ 


1 1 5 4R71Q10 

1 1 J.tO / J ;? 1 VJ 


R91 41Q94 

OZJ .T-JyZT 


9 1 R 

Z. . 1 o 


116.0 


827.25723 


827.25718 


0.29 


116.4873910 


830.88663 


2.28 


117.0 


834.70543 


834.70552 


0.31 


117.4873Q10 


838.33915 


2.18 


118.0 


842.15352 


842.15380 


0.32 


118.4873910 


845.78408 


2.18 


119.0 


849.60183 


849.60199 


0.27 


119.4873910 


853.23131 


2.08 


120.0 


857.05060 


857.05057 


0.35 


120.4873910 


860.68308 


2.18 


121.0 


R64.49914 


R64.49900 


0.30 


121.4873910 


868.13074 


2.08 


122.0 


871.94730 


871.94713 


0.25 


122.4873910 


875.57538 


2.18 


123.0 


879.39555 


879.39559 


0.26 


123.4873910 


883.02764 


2.18 


124.0 


886.84395 


886.84387 


0.27 


124.4873910 


890.47404 


2.08 


125.0 


894.29242 


894.29236 


0.25 


125.4873910 


897.92584 


2.08 
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Table S4: (continued) 



1 96 


Q01 74076 


Q01 74070 


11 

\J. J 1 


196 4871Q10 

lZ.VJ. i ~rO / Jy lu 






127 


909 1 8996 


909 1 8937 


27 

vy . z. / 


127 4871910 

1 Z. / .to / jyiU 


912 82074 

y j_ z..uz. vy / t 


2 28 


128 

J. Z. Cj . Vy 


916 63756 


916 63783 

y i u. vj j / o J 


35 

VyZ J J 


128 4871910 

± z-o ."O / Jyiu 


920 26855 


2 1 8 

Z.. ± O 


1 9Q 


Q94 08561 


Q94 085Q8 

yZ.T.UOJyu 


U.JJ 


log 4R71Q10 

Iz/.tO / Jy lu 


Q97 7141 1 

yjj 1,1 It! 1 


2 1 8 


1 10 








1 10 4R71Q10 

1 JU.tO / JylU 


Q15 16675 

yJJ . 1 VJVJ / J 


2 1 8 

ZJ. _I_ O 


111 


Q18 QR1 11 


Q18 Q8991 


11 

U.JJ 


111 4R71Q10 

1 J 1 . i_ rO / JylU 


Q49 61 985 


9 1 8 

ZJ. _I_ O 


132.0 


946.42970 


946.43047 


0.33 


132.4873910 






133.0 


953.87853 


953.87896 


0.30 


133.4873910 


957.50762 


2.28 


134.0 


961.32731 


961.32765 


0.31 


134.4873910 


964.95949 


2.68 


135.0 


968.77609 


968.77621 


0.33 


135.4873910 


972.40639 


2.18 


136.0 


976.22531 


976.22515 


0.33 


136.4873910 


979.85442 


2.28 


137.0 


983.67397 


983.67372 


0.34 


137.4873910 






138.0 


991.12232 


991.12194 


0.31 


138.4873910 






139.0 


998.57040 


998.56992 


0.33 


139.4873910 


1002.20122 


1.98 


140.0 


1006.01835 


1006.01771 


0.27 


140.4873910 






141.0 


1013.46672 


1013.46612 


0.33 


141.4873910 
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Table S5: Times of planetary transits. 



cycle # 


measured 


uncertainty 


duration 


model 


model duration 


note 




time^. 


(minute) 


(hour) 


time 2 


(hour) 




Planet b 


1 


30 70061 


I 1 Q8 

I I .yo 


4 80 


30 81 466 
- ju.o if do 


A ^0 

f . JU 




9 o 

Z.U 








16 971 49 

lO.Z / If Z 


3 81 


UtiUl gap 


3 o 

J.U 


6^ 94000 


30 00 


6 79 

O. / Z 


6S 94496 

VJ J .Zf f ZD 


6 93 




A 
f .U 


1 1 9 ^3000 


30 00 
JU.UU 


3 60 
J.OU 


1 1 9 ^4^69 
1 1 Z. Jf JOZ 


3 ^3 




J.U 


1 60 04000 


QO 00 


I 1 OO 

I I .uu 


1 60 00608 
lOU.s'UOs'o 


1 90 




6 o 


908 8494S 

ZUO.of Zf J 


S 30 


3 60 


908 84109 

ZUo.of 1UZ 


3 ^0 

J.JU 




7 


9S6 3S001 


60 00 


8 1 6 


9S6 3998S 


7 80 




8 o 

O.U 


30S 1 3831 


70 


3 84 
J. of 


30S 1 9306 

JuJ . 1 ZJ7D 


3 71 
j. / 1 




Q 
y .vj 


3S9 96001 


30 00 

JU.UU 


S 98 


3S9 9S993 

J J Z. Z J ZZ J 


S 00 
J .uu 




10 


401 35165 


8 06 

(J . v./ V.J 


4 56 


401 35574 


4.24 




11.0 


448.43933 


19.12 


4.08 


448.42355 


3.94 




12.0 


497.42072 


60.00 


5.76 


497.46506 


5.50 




13.0 


544.74023 


10.37 


3.60 


544.68713 


3.48 




14.0 


593.25055 


2.30 


9.12 


593.26556 


8.77 




15.0 








640.98407 


3.29 


data gap 


16.0 


688.61578 


26.96 


9.36 


688.65100 


8.65 




17.0 


737.27374 


10.37 


3.60 


737.27942 


3.31 




18.0 


784.40002 


30.00 


5.52 


784.47644 


5.13 




19.0 








833.53937 


3.58 


data gap 


20.0 


880.63666 


12.21 


3.84 


880.61505 


3.73 




21.0 








929.71033 


4.30 


corrupted data 


22.0 


976.86499 


60.00 


4.08 


976.87207 


3.06 




Planet c 


1.0 


246.64867 


5.07 


5.76 


246.64379 


4.02 




2.0 


550.47591 


5.23 


8.16 


550.47833 


6.12 




3.0 


850.99483 


5.30 


6.96 


850.99053 


6.00 




Orphan 


1.0 


977.363 


5.76 


4.15 









BJD-2,455,000 
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Table S6: Model fitting parameters for the photometric-dynamical model. See the text for 
definitions of the terms. The numbers in boldface refer to the parameters shown in Fig. |S231 



Parameter Name 


Best-fit 


50% 


15.8% 


84.2% 


Mass parameters 










0. Mass of Star A, M A (M Q ) 


1.043 


1.049 


-0.055 


+0.054 


1. Mass ratio, Star B, M B /M A 


0.3473 


0.3462 


-0.0063 


+0.0069 


Planet b Orbit (Epoch BJD 2,454,969.216) 










2. Orbital Period, P b (day) 


49.514 


49.532 


-0.027 


+0.040 


3. Eccentricity Parameter, Jei cos(wb) 


-0.094 


0.000 


-0.075 


+0.067 


4. Eccentricity Parameter, yj~e\> sin(wb) 


0.003 


0.098 


-0.067 


+0.042 


5. Time of Barycentric Transit, 










t b (BJD - 2,455,000) 


-31.367 


-31.353 


-0.010 


+0.011 


6. Orbital Inclination, % (deg) 


89.59 


89.70 


-0.16 


+0.50 


7. Relative Nodal Longitude, AS! 6 (deg) 


0.10 


0.23 


-0.21 


+0.58 


Planet c Orbit (Epoch BJD - 2,455,246.6545) 










8. Orbital Period, P c (day) 


303.158 


303.137 


-0.020 


+0.072 


9. Eccentricity Parameter, y / e^cos(cj c ) 


-0.35 


-0.04 


-0.40 


+0.41 


10. Eccentricity Parameter, »/e^sin(u; c ) 


-0.237 


-0.257 


-0.041 


+0.039 


11. Time of Barycentric Transit, 










t c (BJD - 2,455,000) 


246.985 


246.997 


-0.012 


+0.016 


12. Orbital Inclination, i c (deg) 


89.826 


89.825 


-0.010 


+0.009 


13. Relative Nodal Longitude, Af2 c (deg) 


1.06 


0.99 


-0.50 


+0.49 


Stellar Orbit 










14. Orbital Period, P AB (day) 


7.44837695 


7.44837703 


-0.00000021 


+0.00000021 


15. Eccentricity Parameter, ceb cos(oj eb) 


-0.019778 


-0.019797 


-0.000045 


+0.000044 


16. Eccentricity Parameter, ceb sin(ujEB) 


-0.0125 


-0.0112 


-0.0019 


+0.0019 


17. Time of Primary Eclipse, 










t EB (BJD - 2455000) 


-29.306346 


-29.306342 


-0.000018 


+0.000018 


18. Orbital Inclination, i eb (deg) 


89.34 


89.40 


-0.10 


+0.12 


Radius/Light Parameters 










19. Linear Limb Darkening Parameter for 










Star A, u 


0.4151 


0.4137 


-0.0044 


+0.0044 


20. Stellar Flux Ratio, F B /F A (xlOO) 


0.568 


0.579 


-0.017 


+0.017 


21. Density of Star A, p A (g cm -3 ) 


1.163 


1.176 


-0.025 


+0.024 


22. Radius Ratio, Star B, R b /Ra 


0.3636 


0.3671 


-0.0047 


+0.0047 


23. Planetary Radius Ratio, Rb/R A 


0.0283 


0.0289 


-0.0011 


+0.0011 


24. Planetary Radius Ratio, R C /R A 


0.0439 


0.0440 


-0.0018 


+0.0017 


Relative Contamination, 










100 x (F cont /F A ) 










Season 


-2.9 


-0.9 


-2.8 


+2.8 


Season 1 


-1.5 


0.5 


-2.8 


+2.9 


Season 2 


-2.8 


-0.8 


-2.8 


+2.8 


Season 3 


-1.9 


0.1 


-2.8 


+2.9 


Noise Parameter 










Long Cadence Relative Width, ctlc (x 10 5 ) 


62.95 


62.75 


-0.42 


+0.43 


Radial Velocity Parameters 










RV Offset, 7 (km s" 1 ) 


4.67 


4.60 


-0.23 


+0.22 


Zero-level Diff., A7 (km s _1 ) 


-0.12 


-0.01 


-0.31 


+0.30 


RV Jitter, a RV (km s _1 ) 


0.31 


0.43 


-0.12 


+0.19 



Table S7: Derived parameters from the photometric-dynamical model. 



Parameter 


Best-fit 


50% 


15.8% 


84.2% 


Bulk Properties 










Mass of Star A, M A (M ) 


1.043 


1.049 


-0.055 


+0.054 


Mass of Star B, M B (M ) 


0.362 


0.363 


-0.013 


+0.012 


Radius of Star A, R A (R Q ) 


0.964 


0.963 


-0.017 


+0.017 


Radius of Star B, R B (R e ) 


0.3506 


0.3533 


-0.0063 


+0.0060 


Radius of Planet b, R b (i? e ) 


2.98 


3.03 


-0.12 


+0.12 


Radius of Planet c, R c (-R®) 


4.61 


4.62 


-0.20 


+0.20 


Density of Star A, p A (g cm~ 3 ) 


1.163 


1.176 


-0.025 


+0.024 


Density of Star B, p B (g cm -3 ) 


8.41 


8.24 


-0.20 


+0.21 


Gravity of Star A, log g A (cgs) 


4.488 


4.492 


-0.011 


+0.010 


Gravity of Star B, \ogg B (cgs) 


4.9073 


4.9017 


-0.0067 


+0.0067 


Orbital Properties 










Semimajor Axis of Stellar Orbit, a AB (AU) 


0.0836 


0.0838 


-0.0014 


+0.0013 


Semimajor Axis of Planet b, a b (AU) 


0.2956 


0.2962 


-0.0047 


+0.0044 


Semimajor Axis of Planet c, a c (AU) 


0.989 


0.991 


-0.016 


+0.015 


Eccentricity of Stellar Orbit, e AB 


0.0234 


0.0228 


-0.0009 


+0.0010 


Argument of Periapse Stellar Orbit, co AB (Degrees) 


212.3 


209.5 


-4.4 


+4.1 


Mutual Orbital Inclination, I b (deg) 


0.27 


0.43 


-0.24 


+0.66 


Mutual Orbital Inclination, I c (deg) 


1.16 


1.08 


-0.42 


+0.46 


Eccentricities Constraints 










Eccentricity of Planet b Orbit (95% conf.), e b 


< 0.035 








Eccentricity of Planet c Orbit (95% conf.), e c 


< 0.411 
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Table S8: Predicted transit times for planet b. 



Epoch Tq-2,455,000 BJD Impact Parameter Transit Velocity (i^/day) Duration (hr) 
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± 


1 


039 
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1 


.027 ± 
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6. 


.479 
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1. 
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4. 


.824 
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2. 
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.516 
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0. 
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2. 
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4, 


.574 
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2. 


.176 
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0. 
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1 


.006 ± 


0. 


.499 


11 


.336 


± 


1. 


.781 


2. 


.798 


± 


1 


.204 
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0. 


.731 


1 


.003 ± 


0. 


.504 


9. 


.729 
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1. 


,868 


3. 


343 


± 


1 


.556 
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2754. 


.815 


± 


0. 


.941 





.933 ± 


0. 


.477 


9. 


.607 


± 


2. 


.221 


3 


.583 


± 


1 


.858 
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933 
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.977 ± 





.482 


11 


.443 
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1. 


.205 


2. 


.886 
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1 


.109 
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Table S8: (continued) 



\J\J 


28^0 864 ± 1 21 7 


872 ± 4^3 


7 4^3 ± 2 09^ 

f .^UtJ l Zj . WZj 


4 


872 ± 2 1 82 


61 


98Q8 907 4- fl fifi^i 


U.jOI _1_ U.thJU 


1 nqi _|_ n onn 


9 

Zj . 


sin 4- n 8Q7 


6? 


9Q4fi ^(14 4- 1 ^fi 




U.UI 1 —1— ±.0/0 


(J . 


o«q _|_ 000 

0\J<_/ _l_ Zv.ZvZv*J 


63 


2994.499 ± 0.731 


0.874 ± 0.411 


12.231 ± 1.283 


2. 


.896 ± 0.922 


64 


3041.964 ± 1.214 


0.796 ± 0.400 


6.417 ±2.481 


5. 


.578 ±2.385 


65 


3090.777 ±0.898 


0.799 ± 0.370 


11.615 ± 1.982 


3. 


.255 ± 1.367 


66 


3137.912 ±0.997 


0.778 ± 0.366 


9.062 ± 2.402 


4. 


.034 ± 1.866 


67 


3187.001 ± 1.142 


0.713 ±0.332 


10.124 ±2.354 


3. 


.862 ± 1.931 


68 


3234.085 ± 0.863 


0.728 ±0.329 


11.034 ± 1.705 


3. 


.445 ± 1.234 


69 


3283.097 ±1.385 


0.623 ±0.298 


8.157 ±2.238 


4. 


.894 ±2.199 


70 


3330.344 ± 0.823 


0.663 ±0.292 


11.893 ± 1.143 


3. 


.340 ± 0.832 


71 


3378.853 ± 1.505 


0.548 ±0.267 


6.752 ± 2.073 


6. 


.129 ±2.179 


72 


3426.633 ±0.871 


0.579 ±0.265 


12.139 ± 1.403 


3. 


.476 ± 1.016 
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Table S9: Predicted transit times for planet c. 



Epoch 


T -2,455,000 BJD 


Impact Parameter 


Transit Velocity (ii^/day) 


Duration (hr) 


3 


1154.756 ±0.011 


0.430 ± 0.056 


8.131 ±0.133 


4.735 ± 0.260 


4 


1458.197 ±0.264 


0.397 ± 1.258 


3.369 ±0.217 


12.967 ± 1.000 


5 


1758.963 ±0.024 


0.446 ± 0.057 


7.600 ±0.135 


5.086 ±0.278 


6 


2062.831 ±0.033 


0.407 ±0.114 


7.249 ±0.145 


5.953 ± 0.462 


7 


2363.464 ± 0.076 


0.458 ± 0.066 


4.141 ±0.258 


9.823 ±0.783 


8 


2667.055 ± 0.045 


0.436 ±0.116 


8.195 ±0.150 


5.217 ±0.417 


9 


2970.751 ±0.085 


0.462 ±0.164 


5.117 ±0.316 


9.026 ±0.981 


10 


3271.295 ±0.083 


0.462 ±0.118 


6.923 ±0.195 


6.232 ±0.487 


11 


3575.153 ±0.080 


0.461 ±0.172 


7.787 ±0.236 


5.922 ±0.561 



Table SI 0: ELC model parameters. 



Parameter 


Best fit 


e 


0.0306 ±0.0071 


u (deg) 


226 ± 12 


Ra/cl 


0.05322 ± 0.00068 


Rsja 


0.01935 ±0.00029 




0.5958 ± 0.0035 


i (deg) 


89.69 ±0.16 


x A 


0.30 ±0.13 


VA 


0.38 ±0.27 
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